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SUMMARY 


Theoretical studies were made of the effect of spin on the internal ballistics at 


fixed time of rocket motors with axisymmetric grains and the effect of acceleration on 
the burning rate of metallized and nonmetallized composite propellants. In addition, 
the unclassified literature concerning the effects of acceleration on the burning rate of 
metallized and nonmetallized propellants was reviewed. In the internal ballistics 
studies a relatively general internal ballistics at fixed time computer program was 
written and the effects of spin on a single nozzle motor with a CP grain possessing 
inhibited ends were computed. The internal ballistics at fixed time computer program 
was based on the numerical integration along the port of influence coefficient equations 
for one -dimensional vortex flow developed in the study and an adaptation of Mager's 
results for swirling flow in nozzles. The numerical study showed that the major spin 
induced effects in single nozzle motors with CP grains possessing inhibited ends are 
attributable to acceleration induced burning rate changes and nozzle effects due to the 
swirling flow. In the burning rate studies a theoretical model based on particle burning 
was developed for metallized composite propellants, and a theoretical model based on 
the granular diffusion flame combustion model was developed for nonmetallized com- 
posite propellants. The trends predicted by the theoretical metallized propellant model 
relative to variations in burning rate coefficient and pressure exponent with accelera- 
tion and in burning rate with the direction of the acceleration force were in qualitative 
agreement with the experimental data. The theoretical model for nonmetallized com- 
posite propellant also showed qualitative agreement with the general trends in the 
experimental data where comparison was possible. However, a more detailed compari- 
son with the data of Anderson showed important discrepancies. The causes of these 
discrepancies were not resolved. 
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INTRODUCTION 


In all rocket propelled vehicles, or missiles, a high degree of attitude stabiliza- 
tion must be achieved if the actual flight path is to closely approximate the flight path 
desired. At the present time, three methods of attitude stabilization are employed: 
spin stabilization, aerodynamic stabilization, and controlled stabilization. The first 
two methods are passive because neither attitude sensing nor control systems are 
required; the latter is an active system because these components are required. 

It is well known that an appreciable fraction of the cost of a rocket propelled 
vehicle is attached to an attitude control system. Moreover, the inclusion of an attitude 
control system generally decreases the reliability and ruggedness of the vehicle and its 
payload capacity. Therefore, the use of passive attitude stabilization, when the 
generality offered by an active system is not required, is highly desirable from both 
economic and reliability standpoints. 

The method of passive stabilization to be employed depends, to a large degree, 
upon both the characteristics and the flight path of the vehicle. For example, aerody- 
namic stabilization cannot be employed outside the sensible atmosphere because of the 
inability to generate aerodynamic forces of sufficient magnitude. Therefore, consid- 
erable interest resides with spin stabilization, and a number of rocket propelled vehicles 
currently employ this method of attitude stabilization. 

However, numerous static spin tests! have demonstrated that the performance 
characteristics of a Bolid propellant rocket motor depend upon its spin rate. (1-23) 

Figure 1 illustrates the variation of the pres sure -time history for a solid propellant 
rocket motor with a CP grain configuration. The figure shows that the deviations from 
0 rpm performance increase as the spin rate increases and that, in this case, the 
deviations are large; at 400 rpm, the peak pressure is approximately 66 percent greater 
than the 0 rpm value. Moreover, the deviations are not uniform throughout the burning 
period; the ignition phase is essentially unaffected, but the position of maximum pres- 
sure varies with the spin rate. Since high performance rocket motors must of 
necessity incorporate small factors of safety, it is seen that the failure to include spin 
effects in the design of a rocket motor could lead to motor failure. 

The spin effects observed have been attributed to spin induced changes in the 
following: internal gas dynamics, combustion phenomena, metal/metal oxide retention, 
heat transfer, and deformation of the burning surface. Murphy and Wall(^) have 
analyzed the experimental data of Wall^^) and have shown that the latter factor is, at 
most, a secondary effect. Therefore, the remaining factors must control the magnitude 
of spin effects. 

It is important to note that all of these effects are interrelated since (a) the burn- 
ing rate depends upon the orientation of the burning surface, the local pressure, and 
the mass flux relative to the burning surface; (b) the internal gas dynamics depends 


1. Static spin test is employed to denote a test where the vehicle, or rocket motor, is 
essentially fixed in space but rotating about its longitudinal axis of symmetry. 
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upon the geometry of the burning surface, the local burning rate, and the mass discharge 
characteristics of the nozzle(s); (c) the retention of metal/metal oxide depends upon the 
geometry of the burning surface, the flow patterns within the motor, and the local burn- 
ing rate; and (d) the heat transfer depends upon the flow patterns and metal/metal oxide 
retention. Since the hypothesized causes for the spin effects observed in motor tests 
are interrelated and since motor test data usually consist solely of over -all parameters 
(pressure and thrust-time histories), it is difficult if not impossible to determine the 
exact roles of the participating causes from the test data. Therefore, an understanding 
of spin effects must be founded on knowledge of the individual causes and their interplay. 

It is also important to note that spin effects, in the main, lie outside the scope 
of conventional internal ballistics and grain analysis procedures. This precludes the 
analytical design as well as the analytical evaluation of rocket motors for spin applica- 
tions, thereby constraining design and evaluation to trial and error procedures with 
prototype -sized motors--an expensive and time-consuming process. 

In view of this state of affairs an analytical study program that was aimed at 
exploring spin-induced changes in the combustion processes of aluminized and nonalu- 
minized composite propellants and spin-induced changes in the flow patterns in the port 
and nozzle(s) of rocket motors was proposed to the Langley Research Center of the 
National Aeronautics and Space Administration. This proposal led to the present pro- 
gram sponsored by NASA under Contract No. NAS7-406. The program was composed of 
two main phases; 

Phase I - Analytical study of the internal gas dynamics of 
cylindrically perforated grains. 

Phase II - Analytical study of the effects of radial accelera- 
tion upon the combustion mechanism of composite 
solid propellants. 

In Phase I, attention was focused on the effects of Spin upon the flow patterns in the port 
and nozzle(s) of a rocket motor with a cylindrically perforated (CP) grain and the 
attendant effects upon motor performance. 

The tasks to be performed in Phase I were as follows: 

1. Derivation of the equations for flow in the port of a spinning 
rocket motor by integral methods. 

2. Solution of equations derived under Task 1. 

3. Combine solution of port equations with solution of nozzle 
equations , 

4. Parametric study using above equations to determine influence 
of various parameters. 

In Phase II, attention was focused on the combustion phenomenon itself, and the direct 
effect of acceleration upon the propellant burning rate was sought. Two cases were 
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examined: (a) nonaluminized composite propellants and (b) aluminized composite pro- 
pellants. The tasks performed in Phase II were as follows: 

1. Derivation and solution of equations for modified granular 
diffusion flame combustion model. 

2. Derivation and solution of equations for particle burning. 

3. Comparison with experimental data. 

The above tasks have been completed; the remainder of this report describes the 
technical effort in detail. 
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ANALYSIS 


PHASE I - INTERNAL GAS DYNAMICS 


General 


The objective of this phase was to develop an equilibrium ballistics -at -fixed- 
time computer program that would include spin effects. At equilibrium the operating 
point of the motor is established by the equality of the mass rate of gas generation 
through combustion and the mass rate of gas discharged through the nozzle. The 
characteristics of swirling flow through a converging -diverging nozzle had been pre- 
viously studied by Mager. Therefore, the nozzle aspects of the problem could be 

readily solved by adapting Mager's results to the conditions at hand. On the other hand, 
the flow field within the port was a considerably more difficult problem because no 
preliminary work had been performed. An exact analysis of this problem would have 
required, in the least, the numerical solution of two non-linear third order partial 
differential equations. In view of the time allowed and the introductory nature of the 
analysis, an approximate solution was sought instead. The approach taken was to find 
suitable approximations to the velocity profiles and then combine these profiles with 
the once integrated form of the equations of change. This yielded a system of coupled 
ordinary differential equations that could be numerically integrated along the port. 


The discussion that follows concerns swirling flow in the port of an axisymmetric 
grain, swirling flow in converging -diverging nozzles, and the resulting equilibrium 
ballistics -at -fixed-time computer program. 

Swirling Flow in the Port of an Axisymmetric Grain 

Figure 2 shows the system under consideration and illustrates the coordinate 
system, the control surface, and the nomenclature employed. Conventional cylindrical 
coordinates are employed with r denoting the radius and z denoting the axial distance 
from the head end of the motor. The radius of the burning surface is r Q and the radius 
of the throat is r*. The components of velocity in the r, 0 , and z directions are V r , 

, and V z , respectively, and s denotes the distance along the surfac^of the grain. 

The first major question that arises is the importance of viscous effects. 
Experiments with vortex chambers(^) have shown that as long as the radial Reynolds 
number (Re r = p V r r l\x) is greater than 5, viscous effects are negligible. The order 
of magnitude of the radial Reynolds number is readily estimated since 0(PV r ) = 
r b^p* * Thus, the order of magnitude of the radial Reynolds number is 

O (Re r ) = r b P p r 0 /fi (1) 

Typical values for the terms in Eq, (1) are as follows: r^ 0.3 in/sec; Pp 0.064 
lb /in 3 ; r Q 1 in.; and p » 3,8 x 10"^ lb/in-sec. Therefore, the order of magnitude 

of Re r in a typical rocket motor is about 500, Since this value is approximately 100 
times the aforementioned critical value, it may be concluded that, except for the region 


1. Symbols are defined in the Nomenclature. 
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very near to the axis of rotation, viscous effects are negligible. This same conclusion 
was reached by Reference 28 although a different approach was employed. In that case 
solutions for viscous, incompressible vortex flow were examined under conditions 
compatible with those inside a rocket motor. 

Since viscous effects will be small, it was assumed that the flow was inviscid. 

In addition to this assumption it was also assumed that the flow was steady, axisym- 
metric, and adiabatic, and that the fluid was a perfect gas with constant specific heats 
and molecular weight. With these assumptions the equations of change for the volume 
bounded by the control surface illustrated on Figure 2 became: 


Continuity 

$ p v • d A =0 

(2a) 

Energy 

#pHV’dA = 0 

(2b) 

Axial Momentum 

§p VJ'dA + § pdA x = 0 

(2c) 

Axial Moment of Momentum 

§ p t V Q V’d A = 0 

(2d) 


By performing the indicated integrations over the control surface and then differentiating 
those expressions with respect to z, Eqs. (2) were expanded to the following: 

r (z) 
o 

d[ 2 IT Jp^rdr] =dm ( 3a) 

o 

r (z) 
o 

d |^2 7T J pH V z r d r J = H g dm ( 3b) 

o 

r (z) r (z) 

o o 

d 2 7T J prdrJ + dj^2iT Jp V 2 r ^ r J” “ V ! dm =0 (3c) 

o o 

r o (z) 

d [ 2,r j P v z v e r2dr ] = r o v 0s d " 1 (3d) 

o 

To complete the transformation of the equations of change to a set of simultane- 
ous ordinary differential equations, approximate expressions for the radial variation of 
density, pressure, axial velocity, tangential velocity, and total enthalpy must be chosen 
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so that the integrals in Eqs. (3) can be evaluated. Insight into the nature of the functions 
required was gleaned from the inviscid equations of motion. These equations in non- 
dimensional form are as follows: 


a v + 

r 

a v + 

. v + r 

- 

1_.V 

a r + 

T V " " 

Z 3z + 

+ 

r 

a r + 

dr + 

+ v + 3V S 

z — 

a z + 

V + V + 
+ V r V 0 

r 

= 0 




3v + a v + 

— i- + V + 

Sr + 


<3 V - + 

5. 1 = . iE. 

, + J a + 

d z o z 


(4a) 


(4b) 


(4c) 


Interest lies primarily in the regions where spin effects may be appreciable. 
Therefore, the spin rate will also be appreciable and the tangential velocity of the gas 
at the burning surface will be large. Moreover, the analysis will be restricted to 
motors with large length-to -diameter ratios. Therefore, the order of magnitude of 
the fallowing nondimensional terms can be established: 0(z + ) = 1/6 ; 0(r + ) = 1; 

O(V 0 ) = 1; O (V+) = 1; O (V*) = 6 where O (6 ) < 1. If these orders of magnitude are 
employed for the values in Eqs. (4)*, and the terms of lower order neglected, Eqs. 
(4a) and (4c) reduce to 



+ 

r 




I_ a P + “ 0(6) 

+ * + 
p 9z 


(5a) 


(5b) 


while Eq. (4b) is unchanged. However, the order of magnitude of the individual terms 
in Eq. (4b) is 6 . 

Eqs. (5) show two important results: (1) the tangential velocity field controls the 

radial pressure gradient, and (Z) the radial pressure gradient is much larger than the 


O 



= ° (v V 

0(r + ) 


1. 


= 1; etc. 
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axial pressure gradient ^ O | — > Oj — ■— J J . The first result enables compu- 

tion of the radial pressure variation from the tangential velocity field. The second 
result has broader significance since it implies that, to an approximation of order 6 , 
axial pressure gradients are negligible when compared with radial pressure gradients. 
Therefore, to a first approximation the static pressure could be taken as independent 
of axial distance. 

The radius of the burning surface may vary with distance from the head end of 
the motor i~r 0 = f(z) j; therefore, the kinetic energy of the combustion gas may also 
vary from point to point on the burning surface. However, the flame temperature (and 
static enthalpy) will not vary over the burning surface and the preceding analysis has 
indicated that to a first approximation the static pressure can be considered constant 
over the burning surface. The question is, "How important are the kinetic energy 
changes?" This question can be answered, in part, by examining the ratio of the kinetic 
energy of the gas at the burning surface to the flame enthalpy since this ratio is repre- 
sentative of the importance of the kinetic energy to the whole. This ratio is 


r 


H 


- o 


2c T 
P 


2 


f 


( 6 ) 


From structural considerations based on turbine practice O (r Q 0> ) = 1000 ft/sec, and 
based on current solid propellants O(Tf) = 6000°R and O (cp) =0,6 B/lb^-R. There- 
fore, 0(T h ) = 0.006, The pressure change due to this kinetic energy (velocity head) 
can also be determined since for small changes 


= JL. ~ 1 — r 

p y-lt y -1 H 


(7) 


Thus, since 7 ~ 1.2 for solid propellant gases, 0(dp/p) = 0.036. These results show 

that pressure and total enthalpy changes due to the kinetic energy of the combustion 
gases at the burning surface are small. When these results are combined with the 
previous results, it is seen that to a first approximation both the total pressure and the 
total enthalpy are the same for each streamline in the flow field. In other words, the 
flow field is nearly isentropic. This implies that first approximations for the required 
velocity, enthalpy, pressure and density profiles may be obtained from potential flow 
results. 

It is interesting to note that the streamlines with the highest total pressure 
originate near the head end of the motor. However (in the flow field), these streamlines 
lie at the smaller radii where viscous effects (which will reduce their total pressure) 
are most severe. Therefore, real fluid effects tend in some ways to force the physical 
system toward the idealized system assumed, 

— * 

For a potential flow curl V = 0 so that in an axisymmetric flow field 

5 (rV 0 ) 


n 


(8a) 


10 


a v 

i 

d z 


d V 

z 

a t 


(8b) 


a < r V 

a r 


= 0 


(8c) 


Integration of Eq. (8a) shows that 

r V Q = K (r) 

but Eq. (8c) shows that K(r) = constant so that 


(9) 


v e ■ K ' r 


Eq. (8b) in nondimensional form is 


( 10 ) 


a v 


a z 


a v 


a r 


( 11 ) 


Integrating Eq. (11) from the burning surface inward gives 

+ 


; iv* 

v z - 1 


a z 


+ dr + + (1, z + ) 

* f 3V r ' 

From the previous order of magnitude results O 

this result and Eq. (12) \ a z + / 


= 6 


(12) 

Therefore, from 


V + “ V + (1, z + ) + o [ 6 2 (1 - T*) 1 
z z L -i 


(13) 


Since r + is less than or equal to unity and O (V + ) = 1, it is seen that the axial velocity 
is essentially constant across the port. Accordingly, it was assumed that 


V ^ V (r , z) = U (z) 
z z o 


(14) 


The fact that the flow field is isentropic to a first approximation also implies 


that 


Cp 


(15) 


Therefore, the radial pressure and density distribution can be determined by integrating 
Eq. (5a) using Eqs. (10) and (15). This yields 
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p = p o [‘ ' ^~ T ~ M e - 1 1 ] 


y -i 


( 16 ) 


P *P C L 1 


- 1 




(17) 


where M q = K/(r 0 c o ). 

Equations (16) and (17) show that at some finite radius r c density and pressure 
vanish and that for r < r c both density and pressure become negative. Therefore, only 
the flow in the region r c ^ r ^ r Q can have physical significance and the lower limit of 
integration on Eqs. (3) must be changed to r c . The critical radius r c can be determined 
by equating either Eq. (16) or Eq. (17) to zero and solving for r. This yields 


r 

c 


r 2 ~ 

2 


1/2 


M 


0 


1 

2 e j 


(18) 


In the actual problem the axial variations of the variables are being sought. 
Therefore, it will be assumed that the flow field at any axial station is isentropic. This 
is in agreement with the order of magnitude analysis and the physical problem. It will 
further be assumed that the profiles represented by Eqs. (10), (16), and (17) express 
the radial variations of the respective functions and that the axial variations in 
tangential velocity, static pressure, and static density, can be obtained by assuming 
that K, p Q , and p Q are functions of axial distance. In addition it will be assumed that 
the axial velocity can be represented by 

V z (r, z) = U(z) + u f (r,z) (19) 

and that the total enthalpy is a function of axial distance only. The function u 1 (r, z) in 
Eq. (19) is included to account for effects due to the boundary layer at the head end of 
the motor. 


It is convenient to define as a new variable 7] = (r/r Q ) 2 . With this variable the 

common bracketed term of Eqs. (16) and (17) can be written as 

y ~ l ■> r _ v „ 


1 - 2 M 0 [ (^) - 1 


] ] =1-^ i^L 

-I 1 - T] T) 


f (TJ, V c ) 


where 


% = (r c /r Q ) 2 . 


(20) 
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The above results may be summarized by stating the functions to be employed 
to evaluate the integrals in Eqs. (3) 


Vg = K (z)/r 

(21a) 

V = U (z) + u (r, z) 

z 

(21b) 

y 


r 'i 7 " 1 


p = p 0 (z) L f ^> ^ ,j J 

(21c) 

1 


P = P o (z) [ f (v, V c )] 7-1 

(2 Id) 

H = H (z) 

(21e) 


Examination of Eqs. (3) shows that five different integrals must be evaluated 
using Eqs. (21) in order to complete the reduction of the equations of change to ordinary 
differential equations. 


The first integral is 
r 


1_ 
y "i 


2 ff J p V z rdr = A q P q U j [f (V, V c ) ] dt] + 2 n J p u' rd 


(22) 


The two functions on the right hand side of Eq. (22) can be interpreted to be the mass 
flow due to mean flow and the mass flow due to perturbations introduced by the end wall 
boundary layers. With this interpretation 


77 J pu' 


r d r = 0 


(23) 


because a nonzero result would imply net mass transport across the end wall boundary 
layer. Since the flow field has been assumed to be steady, net mass generation within 
the end wall boundary layer is impossible. Thus, Eq. (22) becomes 

r 

o 

2f J P V z rdr = A 0 P 0 (24) 

o 

where the function ip is defined as 

„1 

* = J [£{77 , 77 c ) 


1 

y -i 


d T) 


(25) 
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The second integral is 


2ir\ prdr=A o Pq J [f ( V , V J ] 7 _1 d r, = 


<P 


( 26 ) 


where the function <p is defined as 


The third integral is 


<P = J [ f (T7 » T7 c ] Y _1 d V 

% 

i 


(27) 


2 n I p V z r d r = A o P o u2 I [ £ ^ ] Y _1 d V +4 «u| pu'rd 


r + 
r 


(28) 


+ 2 


TT J p (u')* 


rdr 


Eq. (23) shows that the second term is zero. Therefore, Eq. (28) reduces to 


2 17 f pV 2 rdr = ApU 2 ||)+2fffp (u') 2 r d : 

v Z O O V 


(29) 


The first term on the right hand side of Eq. (29) is the axial momentum due to the mean 
flow while the second term is the axial momentum due to the boundary layer flow. The 
latter momentum is created by viscous effects in the end wall boundary layer. The 
problem now is how does this momentum flow change along the port. Clearly, analytical 
games could be played by choosing different variations, but how could the correct 
variation be established? Accordingly, a variation that seems physically plausible and 
offers analytical advantages will be chosen; therefore, choose* 


r ° 2 

2 it p (u r ) rdr = constant 


The fourth integral is 


(30) 


r 

o 

2 IT J p H V z r dr 
o 


Vo UH J [ £ (T).17 C >] d V 


T 

O 

2 ttH J pu' r dr (31) 
o 


1. Only the change of axial momentum appears in the equations of change; therefore, 
this choice eliminates the end wall effects from the equations. 
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Eq. (23) shows that the second term is zero. Therefore, Eq. (31) becomes 


2 ft [ pHVrdr = A p U H tf) 
J z oo 


The fifth integral is 


y - 1 

2 IT J P v 9 v z r 2 dr = \P Q UK J [f ( TJ , T? c > ] d T) + 2 it K J pu'dr (33) 

oj] o 

*c 

Eq. (23) shows that the second term is zero. Therefore, Eq. (33) reduces to 

r 

o 

2 1 T f p V V r^dr = A p UK ib (34) 

J 0 z oo 

o 

The ordinary differential form desired for the equations of change may now be 
obtained by substituting Eqs. (24), (26), (29), (32), and (34) into Eqs. (3), This yields 

d [ A q P o U 0] = dm (35a) 

d £ p ^ A^ <p J+d^A^p^U 2 p Q d A V* dm z 0 (35 b) 

d[ A o P Q UH l|)] : H g dm (35c) 

d [ A o P o UK *] : r o V 6s drh < 35d > 


Eqs, (35) together with the equation of state 


p ; p Rt 
o *o o 


the supplementary equation 


H : c t | 1 + (M 2 + M 2 ) "1 

p o L 2 z 0 J 
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where 


M z U/c (38) 

z o 

and the equations defining V Q , , and $ give a sufficient number of equations to define 
the dependent variables along the port. 

The nonlinearity of the equations and the general nature of the boundary condi- 
tions preclude an analytical solution. Therefore, a numerical procedure to solve the 
equations was sought. The procedure employed parallels a successful procedure 
currently employed at Thiokol for one -dimensional flow problems. This procedure 
uses a Runge-Kutta numerical integration technique to perform a forward integration 
of the equation for dM^/M^ from the head end to the aft end of the grain. To employ 
this procedure Eqs. (35) must be put into influence coefficient form. Equations (35) 
were first transformed into the dimensionless form 


d A 


d p. 


U 


(I .i) m 2 ^* — 

v tp A ip z U p 

o o 


$ 


2 <p 


d m 


m 


ib 2 

M 

d m 

<p z 

m 

_ 2 vi 

d m 

2 u 

m 


(39a) 


(39b) 


dH dH 


c t c t 

p o p o 


(39c) 


dV< 


0 o 


dV„ 


1 dA 

_1_ o 

2 A 


(39d) 


where 


dH 

c t 
P o 


H - H 


c t 
P ° 


dm 

m 


(39e) 
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and 


dv e . v e»- v 8o 


d m 


V 


6 o 


6 o 


(39f) 


In addition to these equations the following were derived* from the definitions of M , 
Mg , <P , and 


dU 


M 


= 2 


dU 

U 


dt 


(40) 


d M 


M, 


- 2 


d K 
K 


d A 


dt 


(41) 


d <p 

<P 


d M 


0 


<P y 2 

M e 


(42) 


dM 

**■ = H ~ 

^ M 


e 


e 


(43) 


Also, from the equation of state and Eq. (37) 


dp dp 

r o o 


dt 


( 44 ) 


1. The derivation of Eqs. (42) and (43)is presented in Appendix I. 
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and 


d H 
H 



M 


2 


z 



+ M ft ) 


2 V - 1 2 2 

dM — 7- M^ dM fl 

z 2 e e_ + 

2 V - 1 2 2 2 

M z >+ 2 < M . + M e> m 9 


+ 



(45) 


Examination of Eqs. (39a - 45) shows that there are 15 variables and 10 equations* 
Therefore, choose 5 of the variables, that can be specified, to be independent variables; 
the remaining 10 are to be treated as dependent variables. The independent variables 
are dA Q /A 0 , dm/m, -2V , dm/(Um), dH/(cpt Q ), and d^/Vg qI The only dependent 

variables that r^eed tp be expressed in terms of the independent variables are 
dM^/M^ z , dMg /Mg , dt Q /t 0 , and dp 0 /p Q since the other dependent variables can 
be expressed in terms of these variables by algebraic relationships. 

The required algebraic manipulations were performed, and Table I presents 
the influence coefficients for swirling flow while Table II presents the same influence 
coefficients for one -dimensional flow. * Since =0 and ^ = 1 when 

Mg = 0, it is seen that the swirling flow influence coefficients have the correct limiting 
form as M g approaches zero. It is interesting to note that the function C is always 
greater than or equal to one; therefore, a swirling flow "chokes” at an axial Mach 
number that is less than one. 

The differential equations describing the flow in the port of a spinning rocket 
motor were written in terms of the influence coefficients as* 


d M - 

= M 2 
ds 2 



- + G 
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v e o i 


2 if r r G d A 
o b 1 5 o 

m A ds 

o ! 


(46a) 


dM j H - H 

~ - M e G 2 1 + G 2 3 c~t 

L P o 


+ G 
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V 9s " V 9o' 

V 6o J 


2 V r r 
o b 
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dA 

o 

ds 

(46b) 


1. The effect of the axial momentum carried by the injected mass is usually neglected 
in internal ballistics calculations because its effect is small; it was neglecte^l here* 
G— refers to the various influence coefficients in Table I, i,e. = ($ 7 M /<P" *)/ 



INFLUENCE COEFFICIENTS FOR SWIRLING FLOW 
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do 
A o 

H - H V -V„ \ 2fr r, 

,r r s , o 9s0o l ob + 

°45 

dA o] 

d s 

F o 31 33 c t 34 V Q : m 

; p o do, m 

A 

o 

ds ; 




(46c) 


dt 
o 

H - H V_ - V„ ' 27Tr r 

t ■ r - r s + r 0S 0o i ° b * 

G 35 

dA 1 
o , 

d s 

r t u , O __ 

o41 43 ct 44 V- j m 

; ' p o 0 o > 

A 

o 

ds ! 


(46d) 


d r’n 
d s 


2 TTr r_ 
o b 


(46e) 


Distance along the burning surface was employed as the independent variable because 
it is continuous and monotonically increasing. The burning rate was assumed to be 
related to the flow field variables by the equation 


= a 


k 

p n 1 + k (p M 3 
o L 1 o rex J 


(47) 


where ct n, k^, and are constants and the Mach number of the flow relative to the 
rotating surface is 


M 

re 


1 



v e„- (l o n) 




1/2 


/ c 


(48) 


The velocity components required in Eq. (48) were determined from the Mach numbers 
and the static temperature 





1/2 


U 


(49) 
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v eo=K (TRt „»] 


1/2 


(50) 


The initial conditions for Eqs. (46) are 

2 

M = m = 0 
z 



u\ = (r o O) Z /(yRp 


(51a) 
(51b) 
(51c) 
(5 Id) 


The conditions throughout the port can now be determined by forward integration 
of Eqs. (46) from the initial conditions given by Eq. (51). The integration was performed 
numerically with a fourth order Runge-Kutta procedure. 

Swirling Flow in Converging-Diverging Nozzles 

As mentioned previously a theory for isentropic swirling flow in converging- 
diverging nozzles has been formulated by Iviager . The analysis was approximate 
and proceded from a scalar potential function. The results of the analysis showed that 
both the mass flow and the vacuum thrust were swirl dependent. In particular Mager 
showed that 


m /m i = f (a*, y) (52a) 

and 

F/F 1 _ d = g(a *, y, e) (52b) 


where O' * is a swirl parameter given by the expression 


K* n fy -1 

= A "2 


( 53 ) 
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Figure 3 . Schematic Diagram Illustrating Systems for Nozzle Analysis 
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and ^ and Fj_q are the mass flow rate and thrust without swirl. 

The extension of Mager's theory to the problem at hand was accomplished by a 
moment - of-momentum balance between the exit of the port and the throat of the 
nozzle(s). Figure 3 is a schematic diagram that illustrates the systems employed for 
the analysis and shows the nomenclature. It was assumed that in the region between 
the exit of the port and the throat of the. nozzle(s) that the flow was isentropic and that 
moment -of-momentum was conserved (MM)^ = (MM) j J. The first assumption is 
compatible with existing internal ballistic procedures. The second assumption is 
obviously true for axisymmetric geometries where external torques are limited solely 
to those produced by shear stresses; in multiple nozzle configurations this assumption 
is largely unjustified because torques produced by pressure forces are possible. 
However, secondary vortices, which are compatible with the assumption, have been 
observed in test motors. 

The axial moment -of-momentum transported through the control surface at sta- 
tion 1 is 


(MM) 


a 

»j 0 p v 


V r 2 dr 

z 0 


(54) 


The quantities in the integral being evaluated at z = L. The tangential velocity is 
Vg = K/r; therefore, Eq. (54) becomes 

(MM) 1 = Km (55) 

The factor can be computed from the tangential Mach number and the static tempera- 
ture at the port exit (z = L) by employing the definition of the tangential Mach number or 


K 

1 


r o (L) [Mg(L) yR t o (L)] 


1/2 


(56) 


Axial moment -of-momentum is transported through station 2 by two modes; 

(1) rotation of the gas in a nozzle relative to that nozzle and (2) rotation of the gas in a 
nozzle relative to the centerline of the motor by virtue of the rotation of the nozzles. 
The moment -of-momentum transported by the first mode is 


(MM) 1 
2 



( 57 ) 
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The moment -of -momentum transported by the second mode is 


(MM) 1 ' 

2 



Hr,, m. 
c 1 , 1 1 


( 58 ) 


The total moment-of-momentum is the sum of these two quantities or 


(mm >2 = L 

i - 1 


K. 

1 


n 


’ c 1, i 


(59) 


In most practical systems the nozzles are identical so that r , . = r . m. = m. 

and Kj = Kj . Therefore, Eq. (59) reduces to c » 1 cl, j* i j* 


(MM) 


m 



Q + K 


(60) 


The K factor at the nozzle throat can now be related to that at the exit of the port 
by equating Eqs. (55) and (60). This yields 


* 2 
K r K — r Q 
1 c 1 


( 61 ) 


Equation (61) shows that for a single nozzle motor (r , = 0) K, = K* However 
for a multiple nozzle motor K* can be equal to, greater than*; or less than ze’ro 
depending upon the relative magnitudes o| Kj and r^ . Since swirl effects in the 
nozzle are dependent upon the value of a and hence K*. Eq. (61) suggests that swirl 
effects can be eliminated by proper choice of r Q ^ 


By employing Eqs. (53), (56), and (61) the swirl parameter a* may be computed 
thereby permitting Mager's theory to be employed for the determination of swirl effects 
in the nozzle(s) of a spinning rocket motor. 

Equilibrium Ballistics at Fixed Time Computer Program 

The preceding sections contain the theoretical background for determining the 
flow fields in the port and the nozzle(s) of a spinning rocket motor. By employing this 
knowledge and a mass balance, the equilibrium operating point was determined. The 
basic computational procedure was as follows: (1) estimate the head end pressure (p 1 ), 
(2) integrate Eqs. (46) along the port to determine the mass flow, K factor, and 
stagnation conditions at the exit of the port, (3) compute the mass flow rate through 
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the nozzle by employing Eqs. (53) and (61) and Mager's theory, and (4) iterate on p 1 
until the mass flows balance. 

Figure 4 is a flow diagram that illustrates the major elements of the computer 
program. The program is divided into eight major blocks: Block 1 - data input, output, 

and conversion to program units; Block 2 - computation of tables of the ( p 9 
functions with Mq as argument and a table of r Q and z with s as argument; Block 3 - 
estimation of p Q (o); Block 4 - determination of non-zero starting line values for the 
Runge-Kutta integration; Block 5 - the numerical integration of Eqs. (46) to the end of 
the port; Block 6 - computation of the mass rate of flow through the exhaust nozzle(s); 
Block 7 - error detector and iterator; and Block 8 - output of calculated results. 

Block 1 is largely self-explanatory. Units conversion transforms the units 
employed for the input data to the ft, slug, sec, R, ft-lb£, lb£ units employed 
throughout the program. 

In Block 2 tables are formed so that functions may be determined by interpolation. 
The cp, and tables are formed from the functions defined by Eqs. 

(25, 27, 1-9, 1-10); the integrals are evaluated by employing Simpson’s rule. ( 30 ) The 
contour of the burning surface is described by a collection of points specified by their 
corresponding r and z coordinates. An additional INFO parameter is specified for 
each point. This parameter determines where the following segment of surface is 
combustible or not. The collection of points is employed to divide the burning surface 
into segments that are truncated cones and the distance from the head end (s coordinate) 
of each point is computed from 


1/2 


s. 

i 


J i-1 


L (r i " 


r i-i> 


<Z. - z. ,) 
1 1-1 


( 62 ) 


The points defining the surface also define the locations where output data are printed. 

The head end pressure is estimated in Block 3 under the assumptions that pres- 
sure is uniform throughout the port and that there are no total pressure losses, swirl 
effects in the nozzle, or erosive burning. 

Experience with a previous one -dimensional computer program where the 
influence coefficient equations were integrated numerically showed that computational 
problems resulted if the Runge-Kutta procedure was started at s = 0. The difficulties 
were associated with zero values for M 2 Z and m. Accordingly, a small step was made 
under the assumption that p Q = p Q (o) to obtain non-zero values for these terms, and the 
Runge-Kutta integration was begun with these values at s = A . The flow diagram 
illustrates the equations employed to determine the flow variables at s = A , 
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Figure 4. Flow Diagram 
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The Runge-Kutta procedure employed in Block 5 was fourth order and was 
available as a subroutine. The subroutine automatically selected the step size 

based on an accuracy criteria. The integration proceeds from s^ to s The depend- 
ent variables at S£ are stored in a table for output purposes. 

In Block 6 the mass discharge through the nozzle(s) is computed. The total 
pressure and total temperature relative to the nozzle(s) is computed first, C is cor- 
rected for temperature deviations, and the mass flow under one -dimensional conditions 
is computed. The K factor and the nozzle swirl parameter are then computed, and 
Mager's theory is employed to find the correction to the mass floiv due to swirl. 

In Block 7 the absolute value of the relative error in mass generation and mass 
discharge is compared with an acceptable error. If the error is sufficiently small, 
control passes to Block 8. However, if the error is unacceptable the head end pressure 
p 1 is re-estimated. For the first three passes through the re -estimation procedure, 
the next head end pressure is determined from 

pj = pj.i < 1 V™d ) j-l j = 1 - 3 (63) 

For the fourth and successive passes the p^ -error table is examined to determine the 
pressure for which the error is zero. 

In Block 8 the computed results are printed out. Values are presented at each 
of the points employed to define the burning surface. These values are computed from 
the values stored in Block 5. 

PHASE II - EFFECT OF ACCELERATION ON BURNING RATE 


General 

The objective of this phase was to develop analytical burning rate models for 
composite solid propellants that included acceleration effects. Experimental data 
have shown that increased burning rates in metallized propellants are intimately 
connected with metal/metal oxide that is retained on the burning surface by the accelera- 
tion force. Therefore, the rate controlling mechanisms in an acceleration environment 
must be different for metallized and nonmetallized propellants. Accordingly, different 
burning rate models were developed for each case. For nonmetallized propellants the 
granular diffusion flame model of Summerfield^ 32 ' was modified to include acceleration 
effects; for metallized propellants an analytical model originated by Reference 28 was 
modified to include particle burning, and preliminary studies were made with another 
particle burning model. 

Metallized Propellants 

In Reference 28 a theoretical model for the effect of acceleration on the burning 
rate of metallized composite propellants was advanced. The model is based on the 
hypothesis that a fraction of the metal that is evolved at the burning surface during the 
combustion process is retained and burns there. The energy released by this metal 
combustion increases the energy transfer to the decomposing surface thereby increasing 
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the burning rate. It was further hypothesized that the process is steady-state and that 
the metal particles retained on the surface of the propellant burn there until they reach 
a critical size such that viscous drag overcomes the acceleration force and the particle 
is literally blown off the surface. 

The increase in burning rate was related to the amount of energy released at the 
propellant surface by the combustion of metal particles through the steady state energy 
balance 6 7 


- ( r J ] h = r w Q f (64) 

b b a = o v b M M K 1 

where r b is the burning rate, h v is the energy required to heat up and gasify a unit mass 
of propellant, w M is the mass fraction of metal in the propellant, Q M is the energy 
released to the burning surface in the combustion of a unit mass of metal, and the 
function f represents the fraction of metal mass that must be removed from a distribu- 
tion of particles such that all particles larger than the critical size are reduced to the 
critical size. The left-hand side of the equation represents the amount of energy 
required for a rate increase from (r b ) a = Q to r fa , while the right-hand side represents 
the amount of energy released at the surface of the propellant by the combustion of the 
metal particles. 

Rearrangement of Eq. (64) gives the burning rate ratio explicitly as 


r -i - 1 

v<v. = o L 1 --— -] <“> 


Equation (65) shows that the maximum burning rate that can be achieved is dependent 
upon the parameter w M Q M /h y and that this parameter must be less than one for 
physically meaningful results. 

The most difficult task was to determine the function f. This was accomplished 
by assuming that the particles were spherical and that the particle size distribution 
obeyed a log normal distribution. After performing suitable mathematical manipulations, 
the following result was obtained * ’ 


f = \ [erfc Un(V§ Mm )/a] 

erfc [ ln(§^/? Mm )/a + 3<r ] ]. 


ic I* 9 or / 2 

- (s ,,/b w ) e 
M Mm 


( 66 ) 
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The critical particle size was determined from a force balance between viscous 
drag and inertial forces. With the further assumption of Stokes flow this yielded 



M r , 
_ go_b 


P P 

'a'go 



a 



(67) 


Thus, for a particular propellant the radius ratio becomes, in general, 


E /E 

b Mm 


[[r b /(ap)] /[r b /(ap)] c 



( 68 ) 


For a particular propellant in a particular motor at equilibrium conditions Eq. (68) 
reduces to 




(a c /a) 


1/2 


(69) 


Examination of Eqs. (68) and (69) shows that by specifying C7 , and the 

operating conditions, the burning rate increase due to acceleration can be determined. 

The approach presented above in brief represents the first serious and thoughtful 
approach to the problem. Unfortunately, the approach deviates from observed physical 
phenomena. First, the burning surface is pitted. This indicates that increased 

burning rates are attained only at certain locations; the theory assumes the rate increase 
is uniform. Second, particulate residue whose size is large compared with the particle 
size of the metal additive is retained in the motor after firing. This suggests 

that agglomeration occurs on the burning surface; the theory neglects agglomeration. 
These two discrepancies together with the fact that the theory gives no information 
regarding the effect of non-normal acceleration forces form the basic objections to the 
theory developed by Reference 28. 

The basic objections listed above are largely concerned with the neglect of 
particle burning effects. Accordingly the theory of Reference 28 was modified to include 
particle burning effects within the general framework of the theory. In particular, the 
following assumptions were introduced; 

1. Condensed phase particles that are initially retained on the burning 
surface agglomerate and remain on the burning surface. 

2. The interaction between the retained condensed phase material and 
the burning surface occurs only at a finite number of points. 
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3. The line of descent of the agglomerated metal through the propellant 
is colinear with the acceleration force vector. 

The other assumptions were unchanged. 

Figure 5 illustrates the analytical model employed. The figure shows conical 
pits in the burning surface, the mean burning surface, the agglomerated metal particle, 
and illustrates some of the nomenclature. 

The rate of descent of a single agglomerated particle through the propellant 
( r ba) was computed from an energy balance between the amount of energy required to 
increase the burning rate above the base rate and the energy supplied through combus- 
tion of metal at the particle. The energy required to increase the burning rate above 
the base rate is 


E r = [ r ba- (r b ) a=o] h vPp A i < 70 > 


where A { is the effective interaction area of the particle with the propellant. This area 
was assumed to be the area where the local burning rate was modified (refer to Figure 5) 
and is 


A. 

1 


T7 COS 

a 


(71) 


However, co S 2 9 C = 1 -si n 2 9 c and sin 6 c = (r b ) a = Q /r ba . Therefore, Eq. (71) becomes 



(72) 


Since the surface is pitted, any metal added to the agglomerated particle must 
originate within its own cone. If it is assumed that there are N s pits per unit area of 
the mean burning surface 1 , the area swept by each pit is 


A 

s 


- N - 1 CO g § 


(73) 


1. This assumption is equivalent to relating the number of pits to the physical structure 
of the propellant. 
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Therefore, the metal added to the agglomerated particle in unit time is 


r -i p _ 

ba'p M 


A. + (A 


A.) 

1 


( 74 ) 


It has been assumed that all of the metal particles evolved in the interaction area are 
agglomerated but that only the fraction G of the metal evolved on the sides of the cone 
is agglomerated. Since the process has been assumed to be invariant with time, this 
metal must be burned in unit time. However, all of the energy released in this process 
will not reach the interaction area. To compute the fraction that does, it was assumed 
that the energy release was uniform over the surface of the sphere. Therefore, if A 
is the lower heating value for the combustion of the metal in the products of combustion 
of the binder and the oxidizer, the energy supplied to the burning Burface is 


E 

s 


2 77r n § 0 w , 

ba a H p M 


AH 


;a - a.] 

1 s If 


][>-i 


i r bi 



(75) 


The burning rate r^ a may be determined by equating Eqs. (70) and (75). This 

yields 


ba 

! r u 

\ bj a = o 


1 -11(1 - G) 

i - n [l + ip 


1/2 


i) g ]* 


(76) 


2 2-1 

where H z AH^ / (2 h v ) and (3 - (N g 77 C & ) cos Equation (76) gives the 

rate of descent of the agglomerated particle through the propellant. However, the burn- 
ing rate desired is the regression rate of the mean burning surface r^. Examination of 
the vector diagram on Figure 5 shows that 


r - r cos § (77) 

b ~ ba 


Therefore, the burning rate ratio desired becomes 


1 - II (1 - G) 


1/2 


~ b a = o 


i - n [i + ( p - i ) g] 


cos $ 


(78) 
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Equation (78) shows the following: the maximum burning rate occurs when the 

acceleration force is normal to and into the burning surface; the burning rate reaches 
an upper bound at high accelerations (G = 1) given by 



(i - n /3) 


(79) 


so that the term II j3 must be less than one, and some critical angle § c exists 1 
beyond which acceleration does not effect burning rate. The latter occurs because the 
burning rate ratio must be greater than or equal to one. Thus, for angles greater than 
, Eq. (78) is invalid. 

The function G is the fraction of metal evolved on the sides of the cone that is 
agglomerated and retained. Therefore, G is simply the first term of Eq. (66) or 


G 


erfc [ln(S M /? M J > * ] 


(80) 


Moreover, analysis shows that when the drag coefficient is that for Stokes flow the 
radius ratio given by Eqs. (68) and (69) remains valid for particles on the surface of 
the cones. Therefore, the effect of acceleration on burning rate can be determined by 
employing Eqs. (68), (78), and (80). 

Nonmetallized Propellants 

The objective of this part of the study was to extend the granular diffusion flame 
burning rate model developed by Summerfield(^) so that acceleration effects would be 
included* 

The basic features of the granular diffusion flame model are as follows: 

1. A quasi - steady gaseous flame exists adjacent to the burning 
surface . 

Z. The burning surface is dry (the oxidizer and fuel vapors are 
liberated directly from the solid phase by sublimation or 
pyrolysi s) . 

3. No significant chemical reactions occur outside the gaseous 
flame zone adjacent to the burning surface. 


1. $ c is defined as the angle at which r^/f r^) a = 0 


1 . 
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4. The vapors of oxidizer or fuel, or both, are released in the form 
of pockets with a certain mass content and these pockets proceed 

to burn in the surrounding medium of the opposite reactant (granular 
diffusion flame). 

5. Burning occurs as a result of energy fed back from the flame to the 
exposed surface of the propellant, 

6. Transport phenomena in the gas phase reaction zone are molecular 
in nature. 

Figure 6 is a schematic diagram of the granular diffusion flame model. The 
important items illustrated are the heterogeneous solid propellant composed of a solid 
fuel binder and small particles of ammonium perchlorate oxidizer, the burning surface, 
the gas phase reaction zone, and the reacted gases. 

An energy balance at the burning surface was employed to relate the burning rate 
of the propellant to the rate of heat transfer from the gases to the burning surface or 
(considering unit area) 


q 


s 


r b p p 



- T.) - Q I 
l s J 


(81) 


The rate of heat transfer was related to the temperature difference across and 
the thickness of the gas phase reaction zone by 


q" = X ( T - T ) / 6 (82) 

s g f s r 

The thickness of the gas phase reaction zone was related to the nature of the zone 
by solving this problem for limiting cases of low pressure and high pressure and then 
joining these solutions. At low pressure, since the rate of molecular diffusion is very 
much faster than the oxidation reaction, the zone was assumed to be a premixed flame 
in which a second order reaction was taking place. The thickness of the gas phase 
reaction zone in this case was shown to be 


i 


6 

r 



(1 - e ) 2 p A exp (-E/RT ) 1 
L g g J 


(83) 


On the other hand, at very high pressure, the chemical reaction rate is so fast 
that the rate of oxidation is controlled entirely by the rate of inter diffusion and the 
burning rate is determined by the granular nature of the zone. Since the pockets of fuel 
vapor are embedded in the oxidizer vapor, the thickness of the gas phase reaction zone 
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Figure 6. Schematic Diagram Illustrating the Granular Diffusion Flame Model 
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was taken to be the product of the lifetime of an average pocket and the average velocity 
of the gas in the zone or 


6 Z V T r 

r g f v 


(84) 


The average gas velocity was related to the burning rate by a mass balance at the burn- 
ing surface or 


v = t kP /P (85) 

g b'p g ' ' 

The average lifetime of a fuel vapor pocket was related to the characteristic dimensions 
of the pocket and the diffusivity of the gases by 


r fv “ 4/ D g < 86 > 

The mass of the pocket of fuel vapor was related to the mean dimension of the pocket 
and the mean density of the gases by 

m « p g d 8 v (87) 

Therefore, by combining Eqs. (84-87) and rearranging, it was found that 


6 " 


OC 


r p m 


5/3 


2/3 


D 

g 


( 88 ) 


The thickness of the gas phase reaction zone at intermediate pressures was 
determined by a linear combination of the thicknesses for the limiting cases or 


6 

r 


z. 6 + 6 

1 r 2 r 


(89) 


Introduction of the above expressions into Eq. (81) yielded 


(r )'\ 
o a = 


~ P Z [ C (T - T.) - Q 
p p S 1 s 

1 

2 

r 1/3 

1 Z 2 m 

X (T - T ) 
g f s J 


p A exp (- E/R T ) + ,1/2 5/6 


L g u g D p 

g ^ g J 


( 90 ) 
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The scope of this work is the pressure range of practical interest (p > 600 psia). 
In this range, experimental data presented by Reference 32 show that the combustion 
phenomenon is dominated by the rate of interdiffusion; this result is consistent with the 
basic analytical model. Therefore, the analysis will be limited to the effect of accelera- 
tion on the limiting case for high pressures. 

With the granular diffusion flame model, acceleration induced changes in burning 
rate must originate from acceleration induced effects within the gas phase reaction 
zone. Two effects appear possible. First, the acceleration field will produce a pres- 
sure difference across the gas phase reaction zone. Second, the acceleration field will 
act on the density inhomogeneities within the gas phase reaction zone. 

The order of magnitude of the acceleration induced pressure change across the 
gas phase reaction zone is readily estimated. Figure 7 is a sketch that illustrates the 
system for estimating this pressure difference. Application of the momentum theorem 
to the control surface illustrated yields (employing D'Alembert's principle) 


P s - p f =r b p p (v f- v s ) + 6 r p g a (91) 

Equation (91) shows that the pressure difference is due to two factors: the 
momentum change of the gas and the acceleration induced body force. It is also seen 
that the first term is not explicitly dependent upon the acceleration level. 

The order of magnitude of the two terms in Eq, (91) will be estimated for the 
typical conditions tabulated in Table C-Iandthe maximum acceleration level encountered 
in practice (50, 000 g). This yields 


O [r, p (v - v )] = 2.4 x 10 ^psi (92) 

bp i s 

0[6p a] < 8xl0^psi (93) 

r'g 


Since the pressure level of interest is greater than 600 psia, it is seen that both 
of these effects are negligible. Therefore, pressure is essentially constant across the 
gas phase. In addition, it may be concluded that the observed effect of acceleration on 
burning rate must be caused by the action of the acceleration field on the density inhomo- 
geneities in the gas phase reaction zone. 

Density inhomogeneities within the gas phase reaction zone arise from the follow- 
ing main sources: the heterogeneous nature of the zone and the mean temperature 
gradient through the zone. Thus, the gas phase reaction zone may be pictured as having 
mean characteristics that are dependent only upon the distance from the mean burning 
surface and heterogeneous characteristics that are dependent upon both spatial location 
and time. The latter are typified by pockets of fuel vapor embedded in oxidizer vapor. 
The action of the acceleration field on the mean characteristics of the zone will produce 
a mean fluid motion that will be closely related to free convection flows. On the other 
hand, the action of the acceleration field on the heterogeneous characteristics will 
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produce relative motion between the pockets of fuel vapor and the surrounding oxidizer 
vapors. These two effects will be considered separately. 

The nature of the effects produced by an acceleration field acting on the mean 
density gradient within the gas phase reaction zone will depend upon the magnitude and 
direction of the acceleration force relative to the burning surface. When the accelera- 
tion force is normal to and into the burning surface, no effects due to this source can 
occur because the acceleration force is aligned with the density gradient and no flow 
can result. On the other hand, when the acceleration force has a direction opposite 
to the density gradient, free convection flow should be possible. However, a one- 
dimensional flow is impossible because the denser gas would have to ascend at the 
same time that the lighter gas was descending. Experimental and theoretical results 
show that for the case where the fluid is bounded by impermeable walls no flow occurs 
as long as the product of the Grashof^ and Prandtl numbers is small. However, 

when this product reaches a value around 1700, a cellular flow in which the cold, dense 
fluid ascends in the center of each cell and. the hot, light fluid descends along the rim 
of the cell is initiated. ( 33 ) The effect of this flow is to effectively increase the rate of 
heat transfer from the hot to the cold plate. Reference 33 also shows that this increase 
can be accounted for by a fictitious thermal conductivity and that this conductivity is 
a function of the product of the Grashof and Prandtl numbers. 

These results suggest that a similar phenomenon could occur in the gas phase 
reaction zone when the acceleration force vector is away from the burning surface 
and sufficiently large. However, the question arises as to the effect of the vertically 
upward flow that occurs in the case of burning propellant. Since increased heat transfer 
(and increased burning rate) depends, in the case of a cellular flow, upon the downward 
convection of hot, reacted gases, it appears that a vertically upward flow would block 
this flow and delay any appreciable effect of the cellular flow on the burning rate until 
the Grashof -Prandtl product is significantly greater than 1700. 

The Grashof number is 


a p 6 3 (p - p,) . . 

'g r 's 't (94) 

° r 6 “ 2 

H 

' cr 


Therefore, the acceleration required to initiate the cellular flow when there is no mean 
upward flow is 


1700 ix 


cf 




P R 6 r < P s -Pf )Pr 


(95) 


1. The characteristic dimension in the Grashof number is the thickness of the zone 
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Substitution of the typical values found in Table C-I into Eq. (95) shows that a c 3000g. 

In certain applications, acceleration levels of 50, 000 g are encountered. However, it 
is expected that cellular flow effects will not be important in practice because of the 
stabilizing effect of the vertical gas flow and the facts that sustained high acceleration 
occur only in spinning rocket motor and grain configurations that possess propellant 
surfaces normal to and away from the acceleration force cannot be supported at high 
acceleration levels. In other words, structural considerations will, in general, prohibit 
a situation where this effect could be significant. For these reasons, this effect will not 
be considered further. 

The preceding discussion shows that acceleration induced mean flows will not 
occur in situations of practical interest when the acceleration force is normal to the 
burning surface. Therefore, it can be concluded that mean flow effects can occur only 
when a component of the acceleration force is parallel to the burning surface. 

When a component of the acceleration vector is parallel to the burning surface, 
a free convection flow having a boundary layer character will result. Exact treatment 
of this problem is difficult. Therefore, an approximate treatment will be made to 
determine the severity of the effects that could result from these flows. Further work 
on this aspect of the problem will depend on the results of this study. 

In the granular diffusion flame model, the thickness of the gas phase reaction 
zone is controlled by the product of the vertical velocity of the gas and the lifetime of 
a pocket of fuel vapor. The effect of a mean flow parallel to the surface on these 
quantities should not be extremely large. Therefore, the thickness of the gas phase 
reaction zone will remain relatively constant. Since this thickness is small and 
phenomena that can effect the burning rate must occur within this zone, the mean flow 
in the gas phase reaction zone must exhibit a boundary layer character. That is, 
variations of parameters along the burning surface will be small compared with varia- 
tions of the same parameters normal to the surface. Moreover, at large distances 
measured along the burning surface from the edge of the propellant, the variation of 
parameters with distance from the edge will be very small. 

The boundary layer equations for laminar free convection flows of a compressible 
fluid are(^) 


B u B u ^7 f B z' B u 7 

f> + a * + 7 0*17 J 


(96) 


and 


5 (P u) l (P y) o (97) 

Bx By 


1. Very high accelerations can occur in gun boosted rockets while the rocket is in the 
gun. However, the duration of these accelerations is very short. 
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Since variations with distance along the burning surface will be small and the order of 
magnitude of u and v will be essentially the same, the order of magnitude of the term 
v bu/ by will be greater than that of u ^u/b x. Therefore, Eq. (96) can be reduced to 


b u 

P v T“ 

by 


*„ ( p -‘V + T 7 


C“ I7) 


(98) 


Since variations with distance along the burning surface are small, Eq. (97) shows that 


p v =* constant 


(99) 


In Eq. (99) the constant is the value at the burning surface r^p^ 
The viscous shear stress is 


T = M 


b u 
b y 


( 100 ) 


Equation (98) can be rewritten in terms of the viscous shear stress as 


1JL 

b y 



T - a x (p - p f ) 


( 101 ) 


At the edge of the gas phase reaction zone (p - p^) = 0; therefore, for the region 

beyond this zone (y a 6 r ) 


u. = .IJ!e t < 102 ) 

9 y fi 


Equation (102) shows the following: 


(a) 

If T (x, 6 ) > o, 
r 

lim t 

y — » oo 

(b) 

If T (x, 6 ) < 0, 
r 

lim T 


y — » GO 


42 


It can be concluded that u (x, 6 r ) = 0, Therefore, at y = 6 


B u/By = 


B u/B y^ 


= 0 


y ^ 6 


(103) 


The lateral body force on a fluid particle crossing the zone will impart an axial 
momentum to it* It may be concluded that u (x, 6 r ) > 0. 

The velocity profile through the gas phase reaction zone will be approximated 
with the third order polynominal 


2 3 

U = C o + C 1 77 4 C 2 71 + C 3 77 (104) 

where TJ = y/ & r . The four parameters are evaluated by employing Eq. (103), the no 
slip condition at the burning surface, and Eq. (98) evaluated at the burning surface. The 
evaluation yields 
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(p - p ) a 6 R e . 
s f xr 6 

p r Re. + 2 

p b 0 
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[1 - r? +77 2 /3] 


(105) 


where Re = r b*p^ r ^ # Thus, the velocity at the outer edge of the gas phase reaction 
zone is 


u (x, 6 ) 
r 


a x 6 r ( P S - Re 6 
' 3r b°p (Re 6 +2) 


(106) 


The question now arises as to how this mean flow can effect the burning rate. 
Within the framework of the granular diffusion flame model the burning rate is directly 
related to the ratio X / 6 r . Thus, acceleration induced mean flow effects must be 
related to Xg and 6 Two effects appear possible: (1) increase in the effective value 

for Xg through turbulence generated by lateral acceleration forces acting on the density 
inhomogeneities in the reaction zone and (2) decrease in ^ r by the mean flow. The 
first possibility is beyond the scope of this analysis; however, it certainly represents a 
phenomenon that could limit the range of validity of the analysis. The second effect 
requires a mechanism for coupling the acceleration induced mean flow with the lifetime 
of a fuel vapor pocket. Figure 8 illustrates a possible coupling mechanism -velocity 
gradient deformation of the fuel vapor pockets. The figure shows that a pocket embedded 
in a velocity gradient is progressively deformed so that the effective length is increased 
and the effective thickness decreased. Since the lifetime of a pocket is related to its 
smallest characteristic dimension, it is seen that velocity gradient deformation can 
decrease the lifetime of the pocket. Therefore, velocity gradient deformation of the 
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Figure 8# Sketch Illustrating the Deformation of a Spherical Fuel Vapor Pocket by a Velocity Gradient 
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fuel vapor pocket provides a coupling mechanism between acceleration induced mean 
flows and the combustion phenomenon. Note that these flows could be produced by 
pressure gradients (erosive burning effects). 

Figure 9 illustrates the analytical model and nomenclature chosen for the 
analysis of fuel vapor pocket deformation by a velocity gradient. The pocket is idealized 
as a parallelepiped. The figure shows the pocket at some time T after its creation 
(solid lines) and at a later time T + At (dashed lines). The pocket at time T + At 
has been displaced so that its lower leading edge coincides with that of the pocket at 
time t . In the time interval A T the distance S increases so that 


A S 


9u A 

— — A t 

9 y 


(107) 


From Figure 9 it is seen that 


tan (oo + A co) = ( E tan W + H — — 

9 y 

However, the tangent of the sum of the two angles is^O) 


At)/S. 


tan (CO -f- A CO) 


ta n CO -f tan A CO 
1 - tan CO tan A C0 # 


Equating these two expressions, expanding, and simplifying yields 


(108) 


(109) 


tan A to d u 2 d u 

— — = — sec w + tan W ~ At. 


( 110 ) 


Expanding the tangent of A W in a Maclaurin series, substituting this result into Eq. 
(110) and taking the limit of the resulting expression as -At approaches zero yields 1 


d w 
d T 


2 
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CO 


d u 

dy ' 
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1. It is assumed here that u = f(y) and that tan ^ and 3 u/3y are finite. 
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Time is related to the normal velocity of the pocket by 


= v> (ii2) 

d T 


The desired end of the analysis is to determine how a velocity gradient affects the 
lifetime of a fuel vapor pocket traversing it. It is assumed that the pocket is consumed 
according to Piobert's law. Then, if S represents the rate of consumption along a 
normal to the surface of the pocket, the time rate of change of the length of the pocket 
is 1 2 


£_£ = i_^ 

d r cos Jj 


(113) 


Since £ (0) - £ (0) and £ ^ ^ , it is seen that the lifetime of the pocket corresponds 
to the time required for £ to go from £ (0) to 0. ^ 

Combining Eqs. (Ill) and (112) and separating the variables yields 


_d_ 

2 

cos 


JO 

to 


dji 

V 


(114) 


To be in accord with the averaging approach employed by Summerfield in the granular 
diffusion flame model, it is assumed that v = constant. Therefore, Eq. (114) can be 
integrated to yield 


tan to = u /v (115) 

Equation (115) shows that to = 0 at the birth of the pocket and increases as the 
pocket moves away from the surface. 

The effect of this distortion on the lifetime of the pocket can be obtained by 
substituting Eq. (115) into Eq. (113). This yields 


d r 


c os [tan (u/v) ] 
2 S 


d 4 


( 116 ) 


1. It is important to note that d4 /d CO = 0. Therefore, 4 changes only when the 
pocket is consumed, 

2. Note that at any time T 9 £ = £ cos w . Therefore, when 4 = £ = 0. 
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If co " 0, then cos 03 =“ 1-03 U and tan 03 " 03. 


Therefore, Eq. (116) becomes 


d T “ 2 ~ (u/v -l^ d i < 117) 

4 S 


Working with means, we have d T “ Tfv and d£ / ( S) fv) a=0 where(lfv) a _ 0 

is the lifetime of the fuel vapor pocket with no distortion effects. Moreover, u (x, r ) 
u(x, y) and v (x, y) 2 v (x, o). Substitution of these approximations into Eq. (117) 

yields 


T, /( T, ) ^ [ 2 - u (x, 6 )/ v (x, o) ] /2 (118) 

fv fv a = o t 


The acceleration levels of importance to this study are somewhat less than 
1000 g. The velocity u (x, £> r ) can be estimated by employing this acceleration, the 
values tabulated in Table C-I, and Eq. (106). This yields u (x, ° r ) 0.4 ft/sec. A 

typical value for v (x, o) is tabulated in Table C-I as v 8 . Substituting these results into 
Eq. (118) yields T fv /(T,) a=0 “0. 9992. It can be concluded that burning rate changes 
produced by mean flow effects are negligible for the acceleration level of interest m 
this study. It can be further concluded that the effect of acceleration on burning rate is 
connected with the heterogeneous nature of the gas phase reaction zone. 

The effect of an acceleration field on the heterogeneous nature of a gas phase 
reaction zone will be typified by the action of an acceleration field, on a pocket of fuel 
vapor embedded in a steady flow of oxidizer vapor. Figure 10 is a schematic diagram 
that illustrates the analytical model. The schematic diagram illustrates the pocket of 
fuel vapor and the forces involved. 

In the framework of the granular diffusion flame model in the high pressure limit 
the thickness of gas phase reaction zone is 


where v fv is the normal velocity of the fuel vapor pocket and T f v is the lifetime of the 
pocket. Since the burning rate is inversely proportional to the thickness of the gas phase 
reaction zone, the effects of acceleration on burning rate must be explainable by the 
effects of acceleration on V£ v and T £ y . 

The density of the oxidizer and fuel vapors are, in general, different. Therefore, 
in an acceleration field, the pocket of fuel vapor will move relative to the oxidizer vapor. 
This movement will affect v fy . Moreover, the direction of the acceleration will deter- 
mine whether the normal velocity of the pocket of fuel vapor is greater than, equal to, or 
less than the velocity of the flow surrounding it. In addition, relative motion between the 
pocket of fuel vapor and the surrounding oxidizer vapor will always increase the rate of 
interdiffusion, thereby decreasing the lifetime of the pocket. Therefore, it is seen that 
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the action of the acceleration field produces two distinct effects. First, the lifetime of 
the pocket is reduced. Second, the upward velocity of the pocket is altered. It is 
important to note that the first effect is a scalar effect while the second is a vector 
effect. 


Assume now that the density of the fuel vapor is greater than the density of the 
oxidizer vapor. Therefore, when the acceleration force is normal to and into the burning 
surface, the vertical velocity of the pocket is reduced (v£ v < Vg) and the lifetime is 
decreased. On the other hand, when the acceleration force is normal to and away from the 
burning surface, the vertical velocity of the pocket is increased (vf v > v g ), but the life- 
time of the particle is decreased. Finally, when the acceleration is parallel to the 
burning surface, the vertical velocity of the pocket is unchanged (vf v = v g ), but the 
lifetime of the pocket is decreased. It may be concluded that an acceleration force 
into the propellant will cause relatively large increases in burning rate while accelera- 
tions away from the propellant should produce small effects, the direction of which will 
depend upon the relative magnitude of the two effects. Finally, an acceleration parallel 
to the burning surface should cause only moderate to small increases. It is important 
to note that the available data indicate this trend. 

The motion of the pocket is governed by the acceleration indpced body force 
(reversed effective force) A A pa and the drag force p g A Vr$l | Virel/2 |. It will be 
assumed that these two forces are equal and opposite. The volume of the pocket is 

A ^ (120) 


and the frontal area is 


A « d 


£v 


( 121 ) 


Therefore, the relative velocity between the fuel vapor pocket and the oxidizer vapors is 


ad. Ap 
I fv ^ 

rel I p g c d 

L & 


- 1/2 


(122) 


and the normal velocity of the fuel vapor pocket is 


v = v - V , cos $ 
fv g rel 


( 123 ) 
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The mean rate of mass transfer is 


m = h_ A AC 
D v 


( 124 ) 


The mean rate is also 


m = m /t 


fv 


(125) 


Thus, the lifetime of the fuel vapor pocket is 


T r = m 

fv rue 


(126) 


The surface coefficient of mass transfer is usually expressed in terms of the 
Sherwood number <Sh = h D d fv /D g ). Therefore, recognizing that the concentration 
difference is Pf v , Eq, (126) can be rewritten in the following form 


T oc - fv i 

fv — Sh 
g 


(127) 


Comparison of Eq. (127) with Eq. (86) shows that the Sherwood number is a generalizing 
parameter for the lifetime of the pocket. 

The thickness of the gas phase reaction zone is obtained by substituting Eqs. (85) 
(122), (123), and (127) into Eq. (119), This yields 
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(128) 


The burning rate can be related to the thickness of the gas phase reaction zone by 
combining Eqs, (81) and (82). This yields (6 u r = 6^ in the high pressure limit) 


X (T - T ) 

g j 


6 p Lc (T - T ) - Q ] 

T p p S 1 S 


( 129 ) 
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Therefore, 
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( 130 ) 
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When the acceleration approaches zero, the burning rate approaches i|s static value 
(r b ) a= o- Therefore, in Eq. (130), the terms in braces must be ( r b ) a _ Q Pp/ 

(P g C^Sh 0 ) where Sh Q denotes the limiting value for the Sherwood number as the 
acceleration approaches zero. Employing this result and expanding Eq. (130) yields 
with some rearrangement 
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(131) 


Therefore, the burning rate ratio can be obtained by employing the quadratic formula. 
However, before performing this operation, define a Grashof number built upon the 
characteristic dimension of the fuel vapor pocket as 


Gr „ = 


p a Ap d_ 

, fv 


(132) 


and a Reynolds number as 
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It is seen that the term in braces in Eq. (131) is exactly C 2 (Gr d /Re£) 
Therefore, since the burning rate ratio must be positive 
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Dimensional analysis suggests that 


C. = function (Re ,) 
a rel 


( 135 ) 


and 


Sh = function (Re Sc ) (136) 

rel g 

where R e re i = P„ V re }df v / ^ g’ The form of the functions will depend largely upon the 
magnitude of the Reynolds number. 

The relative velocity in the case where the acceleration force is normal to and 
into the burning surface of the propellant ( * = 0) cannot exceed v g . A basic postulate 

of the granular diffusion flame model is that the mass of the fuel pocket in the solid 
will be much less than the mass of a mean oxidizer particle. Moreover, the mass of 
the pocket is invariant with pressure; therefore, 


O (d- ) < d 
fv ox 



and the order of magnitude of the Reynolds number Re 


rel 


will be 


(137) 
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Substitution of the values found in Table C-I into Eq. (138) shows that the Reynolds 
number cannot greatly exceed 20. Therefore, for small accelerations, the flow Bhould 
be dominated by viscous forces. 

For the Stokes flow regime, Redfield and Houghton^ 3 ^ have found that the drag 
coefficients for single bubbles in a liquid obey the relationship 

C d - Re re l" 1 <«9) 


They also report that this relationship is obeyed for liquid -liquid drop systems. The 
mass transfer results of Redfield and Houghton are scattered for the low Reynolds 
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number regime. Therefore, a relationship for solid spheres, which has the correct 
asymptotic form as Re rel approaches zero will be employed 

Sh = Sh D + 0.57 Re^’j Sc®' 35 <140) 


where Sh Q = 2. 

The relative velocity can be obtained by substituting Eq* (139) into Eq. (122) or 


a d fv A P 


rel - 


(141) 


Therefore, the drag coefficient is 


C d = 




(142) 


Substitution of Eqs. (140) and (142) into Eq. (134) yields 
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(143) 


Equation (143) should be valid for an acceleration range where Gr d < 10. For 
accelerations that exceed this criteria, Eq. (143) will fail because the drag coefficient 
becomes a function of the size of the bubble as well as the properties of the fluids. 
Unfortunately, neither empirical correlations nor valid theoretical analyses are available 
for this regime. This precludes the analytical development of a specific equation for the 
burning rate ratio for the regime where Gr^ > 10. 
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RESULTS AND DISCUSSION 


PHASE I - INTERNAL GAS DYNAMICS 


The effect of spin on the head end pressure of a rocket motor with a single nozzle 
and a CP grain configuration was investigated. The rocket motor geometry was charac- 
terized by a grain that was 20 inches long and a single nozzle with d* = 1.5 in. The 
effect of acceleration on burning rate was accounted for by employing the assumed varia- 
tion for O' /{O' ) shown by Figure 11. 

t t a=o 7 6 

Figures 12 and 13 show the effect of spin on head end pressure 1 for the motor 
considered when the port diameter is 1. 8 and 4. 0 in. , respectively. Results are pre- 
sented for the following condition: (1) spin induced changes in burning rate, (2) spin 

induced changes in burning rate and nozzle flow, and (3) spin induced changes in burning 
rate, nozzle flow, and port flow. The figures show that for the lower spin rates the only 
noticeable effect is caused by spin induced changes in burning rate and that for all spin 
rates the effect due to spin induced changes in the port flow is small. In addition, 

Figure 12 shows that head end pressure changes caused by spin induced changes in the 
internal flow field of the motor is small- -even at very large spin rates. However, 

Figure 13 shows that spin induced changes in the nozzle flow are very large at high spin 
rates. 


These results presented suggest that at the acceleration levels of interest to this 
study the primary spin effects in single nozzle motors with CP grains are connected 
with acceleration induced changes in burning rate and spin induced changes in the nozzle 
flow and that the former is the most important. 

PHASE II - EFFECT OF ACCELERATION ON BURNING RATE 


Review of the Literature 


A review of the literature has shown that it can be put into the following general 
categories: literature concerning qualitative spin effects and literature concerning 

quantitative spin effects. The literature in the first category consists mainly of pressure- 
time and thrust -time histories of rocket motors. However, because of the unknown 
magnitude of effects due to internal ballistic changes and non-uniform surface regression 
when the burning rate is dependent upon the orientation of the burning surface, these 
data cannot be reduced to the point where the burning rate can be defined as a function of 
its independent variables. Since this study was concerned with the relationship between 
the burning rate and its independent variables, these data are of little value and will not 
be reviewed herein. This literature comprises References 1 through 23. 

The literature in the second category is composed of literature directed at the 
effects of acceleration on the burning rate of double -base solid propellants, metallized 
composite propellants, and nonmetallized composite propellants. The physical struc- 
ture of double-base propellants and, consequently, their combustion phenomena are 


1, These results were obtained with a special computer program that is limited to 
geometries with r Q = constant. 
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Figure 12. Effect of Spin on Head End Pressure When tt r 



Figure 13. Effect of Spin on Head End Pressure When nr > > A 
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different than those of composite propellants. Therefore, literature that pertains to 
double -base propellants will not be reviewed. As mentioned previously, the combustion 
phenomenon of metallized and nonmetallized composite propellants in an acceleration 
field is different. Therefore, the literature pertaining to these metallized and non- 
metallized propellants will be reviewed separately. 

Metallized Propel lants 

The quantitative literature that pertains to metallized propellants that is unclas- 
sified consists of References 28, 31, and 35-37. These references will be reviewed in 
turn. 


Reference 28 has investigated the effect of acceleration on burning rate of a 
metallized PBAN-AP propellant in the range from 0 to 120 g with a single nozzle motor 
that was spun about its axis of symmetry. The motor had a CP grain with 4. 5 in O. D. 
and 0.6 in. web; the ends were tapered to achieve a neutral area web trace. The 
instantaneous burning rate was computed from the chamber pressure via a mass 
balance and was corrected to the base point pressure by assuming that the pressure 
exponent was unchanged by the acceleration. The propellant variables investigated were 
oxidizer particle size, metal additive particle size, metal content, and additive composi- 
tion. The results of the tests showed that burning rate varied with time during each 
test. In general, the burning rate increased with time and tended toward an asymptotic 
limit. However, in a few tests turning rate decreased with time. In one test the burning 
rate decreased through the firing; in several tests the burning rate first increased and 
then decreased as the web was approached. 

It was also found that acceleration effects were reduced by decreasing the oxidizer 
particle size, decreasing the particle size of the metal additive, and decreasing the 
metal content. In addition, it was found that a refractory metal (tungsten) also produced 
appreciable burning rate increases. Visual examination of the residue showed that no 
agglomeration had occurred and that there was no appreciable chemical reaction. 1 
This suggests that acceleration induced burning rate changes can be produced by a 
mechanism other than metal combustion at the propellant surface. Moreover, when the 
theoretical results were compared with the experimental results, it was found that the 
a c required was far less than the theoretical a c for the additive. This indicated that a 
pre -agglomeration of the metal particles took place before escape from the surface and 
final agglomeration in the pit. Figure 14 presents the effect of acceleration on the mean 
burning rate of an aluminized propellant with PBAN binder. 

Reference 31 has investigated the effect of acceleration in the range from 0 to 
300 g on the burning of an aluminized PBAA-AP composite propellant in a special slab 
burner. The slab burner had a web thickness of 0. 5 inch and was mounted on an arm 
centrifuge so that the angle between the acceleration force and the burning surface could 
be varied. The study showed that acceleration affects both the burning rate coefficient 
and the pressure exponent (Figure 15) and that burning rate at any acceleration level 
depends upon the orientation of the burning surface relative to the acceleration forces. 

In addition, it was found that the amount of residue retained in the motor after firing was 


1. Personal communication with P. G. Willoughby. 
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related to the acceleration level. A sieve analysis showed that the mean size of the 
residue was also related to the acceleration level and that the mean residue size 
increased as the acceleration level increased (Figure 16). Moreover, the normalized 
size distribution of the residue is essentially independent of the acceleration level and 
is that for the aluminum additive (Figure 17). 

Reference 37 has investigated the effect of acceleration on burning rate with a 
strand burner mounted on the arm of a centrifuge. The study investigated several 
different aluminized AP propellants with PBAN and Polyurethane binders at accelera- 
tion levels from 0 to 2000 g; the acceleration force was always normal to and into the 
burning surface. Figure 18 illustrates typical results. The results of the tests showed 
that the burning rate increased rapidly as the magnitude of the acceleration increased 
tending to an asymptotic value at the higher acceleration levels. In addition, it was 
found that the pressure exponent of the PBAN propellant was relatively insensitive to 
acceleration* while that of the Polyurethane propellant was quite sensitive. Residue 
was found in the case that enclosed the strand. At acceleration levels above 300 g, the 
residue was present as a single piece whose thickness increased with the acceleration 
level. X-ray diffraction analysis showed that the residue was aluminum oxide. The 
burning rate data for the PBAN binder with 17% of 7, lJLt aluminum was correlated by 
the equation 


r b /(r b ) a=Q = 1 + 0.028 a 0 - 30 (144) 

while the burning rate data for the Polyurethane binder with 17. 7% of 3 1 JLI aluminum 
was correlated with the equation 


V( r b>a=o 


r 

1 + C 1 ap/(r b ) 



(145) 


where C j and are pressure dependent. 

Reference 35 reported test results for small rocket motors with CP and six 
pointed star grains. In the tests the pressure was approximately 600 psi, the accelera- 
tion ranged from 0 to 100 g, and the acceleration force was parallel to the burning 
surface. The test results showed no acceleration effects. 

Reference 36 has reported test results for both small and moderate sized 
rocket motors containing aluminized composite propellant with CP grains. The tests 
were performed on an arm centrifuge and the acceleration vector was either parallel 
to or normal to the axis of the motor. The tests covered the pressure range from 
1, 000 to 2, 000 psi and acceleration levels from 0 to 200 g. Regression analysis of the 
test results showed no effects that were attributable to acceleration. 

Reference 38 presents results derived from micromotor tests on an arm centri- 
fuge with UTX-2649 propellant containing 16% aluminum. The acceleration range from 


1. This supports the data reduction procedure of Reference 28. 
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Figure 18. Effect of Acceleration on the Burning Rate of a Metallized 
Composite Propellant. (Reference 37) 



65 


0 to 800 g was studied at pressure levels of 1, 000 and 2, 000 psi, in all cases the 
acceleration force was parallel to the burning surface. The tests showed no consistent 
change in burning rate with acceleration level. 

The experimental results suggest the following general trends: 

1. Acceleration induced burning rate changes are largest when the 
acceleration force is normal to and into the burning surface, 

2. Acceleration induced burning rate changes are negligible when 
the acceleration force is parallel to or away from the burning 
surface . 

3. Acceleration induced burning rate changes can be reduced by 
reducing the metal content, metal particle size, and the oxidizer 
particle size. 

4. Acceleration can affect the pressure dependency of burning rate 
and this effect depends upon the binder. 

5. Acceleration induced burning rate changes occur irrespective of 
the reactivity of metal additive. 

Nonmetallized Propellants 

The quantitative literature relating to the effects of acceleration on the burning 
rate of nonmetallized composite propellants consists of References 24, 25, and 28 and 
References 38 through 42. References 25, 40, 41, and 42 are classified and therefore 
will not be reviewed herein. However, Reference 25 has been reviewed by Reference 
24; Reference 42 has been reviewed by Reference 31; and the pertinent data in 
References 40 and 41 are for pressures below the range considered here. The remain- 
ing references are reviewed in the following paragraphs. 

Reference 38 has tested micromotors with CP grains on an extended arm centri- 
fuge. In all cases, the acceleration force vector was parallel to the surface of the 
propellant and in the direction of gas flow. The tests were conducted at nominal 
chamber pressures of 1, 000 and 2, 000 psia and at acceleration levels up to 800 g at 
200 g intervals. Two different propellants, UTX-1724 and UTX -2649, were tested. The 
test data were somewhat scattered; however, there was no indication of an appreciable 
effect of acceleration on burning rate. The essentially null effect observed at two 
different pressures indicates that neither the burning rate coefficient ^ nor the pres- 
sure exponent n in the burning rate law r^ = O' p n are strongly dependent on accelera- 
tion in the range of the test. 


Reference 28 has static spin tested small rocket motors with a PBAN propellant. 
The grain configuration used was an axisymmetric internal burner with inhibited wedge 
shaped ends to provide a neutral trace and maintain the burning surface so that the 
acceleration vector is always normal to and into the burning surface. The nominal 
chamber pressure was 600 psia and the acceleration level ranged from 0 to approxi- 
mately 120 g. Figure 14 presents the experimental results and shows that there were 
no appreciable effects. 
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Reference 24 has reviewed the data of Wall< 25 ) which showed that the pressure 
exponent n of the burning rate law r b =a t p n was not influenced by spin, but that the 
burning rate coefficient O', was dependent upon the configuration of the grain for a 
spinning rocket motor. Reference 24 has also correlated the variation of with 
grain configuration found by Wall by employing the following variables: dimensionless 
burning rate coefficient 0^/(0^) © and the^grain complexity factor X = C/(4^ A c ). 

The correlation showed that &/(a^) a=0 = a X where a and b are constants and b is 
negative. Thus, increasing the grain complexity factor reduces mean acceleration 
effects. In addition to the above. Reference 24 showed that the effect of spin induced 
grain deformations on the performance of the motor was, at most, a secondary effect. 

Reference 37 has investigated the effects of acceleration on the burning rate of 
propellants with polyurethane and PBAN binders. The investigation was performed 
with strands mounted in a bomb. The bomb was attached to the end of a 3 -foot long 
centrifuge arm. In all tests, the acceleration was normal and into the burning surface. 
The tests were conducted at pressures of 500, 1, 000, and 1, 500 psia and at accelera- 
tion levels up to 2, 000 g. Figures 19 and 20 illustrate the effects of acceleration on the 
burning rate ratio r b/( r b^a = o the polyurethane propellant at pressures of 1,000 and 
1, 500 psia. The results suggest that there is a very large effect at relatively small 
accelerations (a < 200 g). However, the large effects occurred in the tests where the 

wire employed to ignite the stand broke during the test. The large increases in burning 
rate measured with the wire broken were attributed to the retention of a portion of the 
broken wire on the burning surface. This explanation appears reasonable because the 
burning rate of a metallized composite propellant is more sensitive to acceleration 
than a nonmetallized propellant (Figure 14) and metal and/or metal oxide is certainly 
retained on the burning surface of a metallized propellant. The large deviation between 
the data points where the wire was intact should also be noted. This suggests that con- 
siderable experimental error existed. 

Figures 21, 22, and 23 illustrate the effects of acceleration on the burning rate 
ratio of PBAN propellant at pressures of 500, 1000, and 1500 psia. These results show 
the following: appreciable acceleration effects occur for accelerations as low as 500 g, 
the burning rate is essentially unaffected by accelerations under 100 g, burning rate is 
highly dependent upon acceleration for 100 < a ^ 300 g, and the burning rate appears 
to approach an upper limit at very large accelerations. The insensitivity of burning 
rate to acceleration at accelerations below 100 g is supported by Reference 28. 

Figure 24 is a composite plot of all of the data at 500, 1, 000, and 1, 500 psia for 
which the ignition wire remained intact. The figure shows that the variation of the 
burning rate ratio with acceleration is essentially the same at all three pressures. 

This indicates that the pressure exponent n is not strongly dependent upon the magnitude 
of the acceleration. The insensitivity of the pressure exponent to acceleration level is 
supported by the data of Reference 25. 

In summary, the review of the literature that is pertinent to the effects of 
acceleration on the burning rate of nonmetallized composite propellants has shown 
that no previous theoretical analyses exist and that the available experimental data do 
not adequately cover the independent variables involved and their ranges of interest. 
Moreover, the majority of the data is derived from strand tests (Reference 37). It is 
well known that strand burning rates differ from motor burning rates. Therefore, 
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Figure 24. Comparison of Burning Rate Data at 500, 1000, and 1500 psia (Data Reference 37) 
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motor data confirming the trends shown by the data of Reference 37 are necessary 
before complete confidence can be placed in them. However, the experimental data do 
suggest the following trends. 

1. The effect of acceleration on burning rate is small for accelera- 
tions less than 100 g, 

2. Acceleration effects depend upon the direction of the accelera- 
tion vector relative to the burning surface. 

3. When the acceleration vector is parallel to the burning surface, 
the effect of acceleration on burning rate is small. 

4. When the acceleration vector is normal to and into the burning 
surface, appreciable increases in burning rate can occur. 

5. When the acceleration vector is normal to and into the burning 
surface, the burning rate approaches an upper limit as the 
acceleration level increases, 

6. The burning rate ratio is essentially independent of pressure. 

Comparison of Theory and Experiment 

Metallized Propellants 

The general characteristics of the analytical model will be examined first. The 
theory developed in the ANALYSIS shows that the burning rate ratio is a function of the 
following variables: j8, a, $ and (ap/r b )/(ap/r b ) c . The latter variable 

contains the acceleration term so it will be treated as the independent variable and the 
other variables will be treated as parameters. Figure 25 illustrates the variation of 
the burning rate ratio with (a p/r^) / (a p/r^) c for the parameter values specified. The 
figure shows that the burning rate ratio increases asymptotically to a maximum value 
and that the way the asymptotic value is reached is dependent upon the variance of the 
particle distribution ® . When O is small the asymptote is reached in a steep step; 
when O’ is large the asymptote is reached in a more gradual step. Figure 25 also 
compares the theory developed herein with the theory developed by Reference 28. The 
comparison shows that, for the conditions shown, the two theories are essentially 
equivalent for small values of the independent . variable when CT = 0,46; however, the 
new theory reaches the asymptotic value sooner. Moreover, the theoretical results of 
Reference 28 are not as strongly dependent upon the variance of the particle distribution 
as the new results are. 

Most of the available experimental results present the burning rate ratio as a 
function of acceleration. Therefore, the independent variable a/a c is desired rather 
than the variable (a p/r^) /(a p/r^) c . This transformation can readily be achieved for 

r mstant pressure conditions (p = p_) because under thi« constraint a/a c = 
(ap/r b )/(ap/r b ) c J [_ r b /(r b ) a = 0 J / T (r b ) c /(r b ) a _ „ J . The burning rate ratio at 
critical conditions is determined from the general results (refer to Figure^ 
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Figure 26 illustrates the effect of the parameters II and j3 on the burning rate 
ratio at constant pressure conditions , The figure shows that the general trends in 
the variation of burning rate ratio with acceleration ratio are unchanged and that both 
parameters affect the maximum value of the burning rate ratio. It is also seen that for 
equivalent asymptotic burning rate ratios (equivalent D fi products) the characteristics 
depend upon the values of the individual parameters and that the magnitude of the 
characteristics change depends upon the magnitude of the II 0 product. When the 
Hr product is small, the characteristics are not strongly dependent upon the 
individual ( II or 0 ) values. 

The theoretical results also show that the burning rate ratio is pressure 
dependent. The pressure dependence is conventionally expressed in the form 


r b = C o ( P / P C ) 


(146) 


The effect of acceleration on the constants and Cj will now be determined when 


P 


The theory shows that burning rate can be written as 


— P c * / 

u ) = “Jp/pJ 11 J 

b a-o t c -J • 


r b = a t (p/p c ,Ilf [ (ap/r b >/(ap/r b ) c ] 


(147) 


where the function f is illustrated graphically by Figure 25. Comparison of Eqs. (146) 
and (147) shows that the burning rate coefficient is C 0 = a t f • The burning rate 
exponent may be determined by differentiating Eqs. (146) and (147) with respect 
to p, equating the resulting expressions for r^ 6 r^/5p , and assuming that C 0 is not 
strongly dependent on pressure. This yields 


n + V [( ap/r b )/(ap/r b ) c ]/f 

TTrT (ap/r b )/( a p / r^F/f 


(148) 


where the prime denotes differentiation with respect to the argument of the function. 
Equation (148) shows that the burning rate coefficient has the same characteristics as 
the burning rate ratio. However, r since f, f 1 , and (ap/r^/lap/r^ are all 
greater than or equal to zero and f 1 (o) = f 1 ? 00 ) = 0, the C\ variation will first increase 
from n and then decrease to n as the acceleration increases. 

Figure 27 shows the theoretical variation of C Q and with acceleration when 
p = p c . Comparison of these results with the experimental results of Northam 
(refer to Figure 15) shows that the theory predicts the general trends shown by the 
experimental data. An exact comparison cannot be made between the figures because 
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Northam's data also show that the parameter fi is dependent upon the acceleration 
(refer to Figure 16); the theoretical results shown assume that P = constant. Further 
investigation of the theory when ft ^ constant is required before a quantitative com- 
parison can be made. It should also be noted that the data of Reference 37 show that 
the pressure exponent depends upon both acceleration and binder composition. The 
variation observed with polyurethane binder agrees qualitatively with the theory; no 
significant variation was observed with PBAN binder. Further work is required to 
investigate the dependency of the pressure exponent with binder composition. It 
should also be noted that the theoretical variation of the burning rate coefficient C Q is 
in qualitative agreement with the data of References 28, 31, and 37. 

The general theory also shows that the burning rate varies with the direction of 
the acceleration force relative to the burning surface. This variation was investigated 
for conditions where the acceleration was sufficiently large for the burning rate ratio 
to be independent of the acceleration (saturated conditions). Figure 28 shows the varia- 
tion of the burning rate ratio with angle for cases where [ ri_/(r K )_ ] are 1.5 

and 2. The figure shows that the burning rate ratio decreases smoothly to unity as $ 
increases. A simple cosine variation is included for comparison purposes. The 
figure shows that the theoretical change is more rapid than that predicted by the cosine 
function. The figure also shows that there is a critical angle beyond which, theoreti- 
cally, the burning rate is independent of acceleration. Figure 29 illustrates the varia- 
tion of the critical angle with the maximum value of the burning rate ratio. The figure 
shows that acceleration effects will not exist for $ > 55°. The data of Northam ^1) 

at 100 g show that when = 30°, r^/(r^) a = 0 = 1, Examination of Northam's data 
shows that at 100 g the burning rate ratio is approximately 1. 13. The theoretical value 
of $ c for this burning rate ratio is approximately 26° (refer to Figure 29). Therefore, 
the experimental and theoretical results again exhibit qualitative agreement. 

Nonmetallized Propellants 

The extension of the granular diffusion flame model presented in the ANALYSIS 
has shown that pressure differences and mean flows induced by accelerations are, at 
most, second order effects and that major effects must originate from the heterogene- 
ous nature of the gas phase reaction zone. Moreover, the analytical formulation of the 
granular diffusion flame model has been extended and an equation relating the ratio of 
the burning rate in an acceleration field to the static burning rate to the drag coefficient 
of a pocket of fuel vapor, the direction of the acceleration force vector relative to the 
burning surface, and suitable Grashof, Reynolds, and Sherwood numbers has been 
derived. Because the first results are not amenable to comparison with experimental 
results, they will not be discussed further and the discussion will center on the analyti- 
cal relationships that have been derived. In the first part of the discussion, the 
implications that arise as a consequence of these equations will be presented. In the 
second part of the discussion, the trends shown by the experimental data will be com- 
pared with the trends predicted by the analytical results. 

Equations (134) and (143) have many similar characteristics. Therefore, these 
similar characteristics will be discussed first. 

The fact that the Grashof number appears in the numerator of both equations 
implies that the burning rate ratio can increase without bound as the acceleration level 
increases. However, the experimental data of Anderson^*?) (Figures 21 through 23) 



Theory 
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and simple logic indicate that when ^ = 0, the burning rate ratio approaches an upper 
bound as the acceleration level increases. When the acceleration is normal to and into 
the burning surface ( $ = 0), the acceleration causes the velocity of the fuel vapor 
pocket to lag the velocity of the oxidizer vapors (vf v < Vg). Therefore, when the 
acceleration is sufficiently large, the lag will have increased to the point where fuel 
vapor cannot escape from the burning surface in an unreacted state. Wherefore, the 
gas phase reaction zone should take on a stratified character. In this event, it would 
be expected that further increases in the acceleration level would have little effect on 
the burning rate. 

The discussion above shows that the granular diffusion flame model must fail at 
high acceleration levels because the gas phase reaction zone loses its granular charac- 
teristics. The parameter that should govern the "granularity 11 of the gas phase reaction 
zone is the ratio of the thickness of the gas phase reaction zone to the characteristic 
dimension of the fuel vapor pocket. Therefore, when this parameter is large, the 
individual pockets will be identifiable and the model should apply. However, when this 
parameter is small, the model will not apply. 

The fact that the Grashof number appears in the numerator also indicates that 
reducing the characteristic dimension of the fuel vapor pockets will reduce the effect of 
acceleration on burning rate. The characteristic dimension of the fuel vapor pocket 
has been assumed by Summerfield to be related to the mean diameter of the oxidizer 
particles. Therefore, reducing the mean diameter of the oxidizer particles should 

reduce spin effects. 

The fact that the static Reynolds number appears in the denominator shows that 
anything that increases the static burning rate will generally reduce the severity of 
acceleration effects. Therefore, the addition of burning rate catalysts as well as an 
increase in the initial temperature of the solid propellant should reduce acceleration 
effects. Moreover, since the static burning rate is related to the mean diameter of 
the oxidizer particles in such a manner that decreasing the mean diameter increases 
the static burning rate, it is seen that reducing the mean diameter of the oxidizer 
particles should be an effective way to reduce acceleration effects. 

It should also be noted that the equations predict different burning rates for 
<j> = 0, 90, and 180°. Moreover, the burning rate for 3> = 0 will be the greatest 
followed by that at $ = 90. The physical reasons for this have been discussed pre- 
viously in the ANALYSIS. 

It is interesting to note that the burning rate at $ = 90° is not dependent upon 
Re 0 and therefore, (r^) a = Q . 

The effect of acceleration on the pressure exponent will be examined next. 
However, the dependency of the Grashof and Reynolds numbers on pressure must be 
determined first, ^ The mass of a fuel vapor pocket is given by Eq. (87). Therefore, 
the Grashof and Reynolds numbers can be rewritten as 


Gr . 


a A p m 

2 

U 

r ct 


1. The Schmidt number Scg is not dependent upon pressure. 


(149) 
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and 


Re 


N 


p m 

a = o 2 _ 


1/3 


u p 

g g 


1/3 


(150) 


Now both Ap and p g are directly proportional to pressure and the static burning 
rate at high pressure is proportional to pl/^. However, the mass of the fuel vapor 
pocket, the acceleration, the dynamic viscosity, and the propellant density are inde- 
pendent of pressure. Therefore, the Grashof number Gr^ is proportional to pressure, 
but the static Reynolds number is independent of pressure. Thus, Eq. (143) becomes 



a=o 


p cos $ + 


[(C & p cos 4>) 2 + C ? p 1/2 + i ] 1/2 


(151) 


Since the static burning rate is proportional to p*/3 f it is seen that the pressure 
exponent predicted by Eq. (151) is greater than unity (n ~ 4/3). The exception occurs 

when £> = 90° (n ~ 7/12). 

This result implies that nonmetallized composite propellant could not be fired 
in a spinning rocket motor. Since the latter has occurred numerous times, a basic 
flaw must exist in Eq. (143). Examination of Eq. (134) shows that pressure dependency 
is largely controlled by the ratio Gr^/C^. Thus, the flaw could occur in connection 
with the relation assumed between C d and Re rei . Re-examination of the data of 
Redfield and Houghton shows that for Re re i > 10, the drag coefficient becomes a 

function of the physical properties of the fluids as well as the size of the bubble. More- 
over, in the analysis no specific account was made of the hot, reacted gases that must 
enclose the pocket of fuel vapor. Since the hot, reacted gases will flow counter to the 
direction of motion of the fuel vapor pocket relative to the oxidizer vapors, the relative 
velocity between the pocket and its immediate surroundings will be greater than the 
velocity of the pocket relative to the oxidizer vapors. This would increase the effective 
R e rel* Furthermore, the system at hand is a gas -gas drop system while the empirical 
data being employed apply to liquid -liquid drop and liquid-gas drop systems. The lack 
of the stabilizing effect of surface tension combined with the small viscosities that exist 
in a gas system could greatly alter the drag characteristics as well as the mass trans- 
fer characterisitcs. Therefore, it appears plausible that the fault noted is connected 
with these factors. 

The trends predicted by the theoretical results can be summarized as follows: 

1. Burning rate in an acceleration field depends upon the orientation 
of the acceleration vector relative to the burning surface. The 
greatest change occurs when $ = 0°. When $ = 90° or 180°, 
changes in burning rate should be small compared with the change 
when $ = 0°. 
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2. The effect of acceleration on the burning rate ratio can be 

reduced when $ # 90° by the following: decreasing the 

particle size of the oxidizer, adding burning rate catalysts, 
and increasing the initial temperature of the propellant. 

3. When the acceleration force is parallel to the burning surface 

( 4> = 90°), the burning rate ratio is not dependent upon either 
the catalyst content or the initial temperature of the propellant. 

In addition to the above trends which were derived from Eqs. (134) and (143), a 
logical extension of the basic model to a situation where the acceleration vector is 
normal to and into the burning surface ( $ = 0) and large shows that the flame loses its 
granular character and that there is an upper bound to acceleration induced changes in 
burning rate. 

The trends shown by the experimental data have been summarized previously in 
the REVIEW OF THE LITERATURE. Comparison of the theoretical and experimental 
trends shows that the theory predicts the correct trend for the variation of burning rate 
with respect to the direction of the acceleration vector relative to the burning surface. 
In addition, the basic model predicts the observed upper bound in burning rate when the 
acceleration vector is normal to and into the burning surface. These results tend to 
confirm the analytical model. 

The discussion above pertains to the general character of the analytical results 
and their agreement with the general trends shown by the available experimental 
results. However, Reference 37 has employed the theory to correlate acceleration- 
burning rate data under constant pressure conditions. Figure 30 presents the results 
of the comparison. The figure shows that the general trend of the data is well repre- 
sented with the following important exception: (1) the theory predicts an immediate 

appreciable increase in burning rate while the data show that no increase in burning 
rate occurs until a > 100 g, and (2) a negative Ap is required to fit the data. 

Reference 37 has suggested that the first defect could be overcome by employing the 
concept of a critical acceleration. That is, for accelerations below the critical value 
r b/( r b) a = o = 1 and ^ or accelerations above the critical the theory is employed. The 
second consequence implies the largest acceleration effect will occur with this pro- 
pellant when the acceleration force is away from the burning surface. Experiments 
with metallized composite propellants have shown that no appreciable change in burning 
rate occurs when the acceleration force vector is away from the surface^ 1 ). Since 
particle burning effects cannot occur when the acceleration force is away from the 
surface, these results suggest that the corresponding effect on a nonmetallized propel- 
lant would also be negligible. Therefore, the negative Ap required by the correlation 
is probably not in agreement with experiment. However, additional data are required 
for positive confirmation. 
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Figure 30. Comparison of Theory and Experiment for Nonmetallized Composite Propellant 

(Reference 37) 
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CONCLUSIONS 


The following conclusions may be drawn from the results of the present investi- 
gation. 

PHASE I - INTERNAL GAS DYNAMICS 


1. The major internal ballistics effects in single nozzle motors with 
CP grains possessing inhibited ends that are due to spin are con- 
nected with acceleration induced changes in burning rate and swirl 
effects in the nozzle. 

PHASE II - EFFECT OF ACCELERATION ON BURNING RATE 


Metallized Propellant 

1. Acceleration induced burning rate changes can occur even though 
the metal additive is non-reactive. 

2. Acceleration affects both the burning rate coefficient and the 
pressure exponent. 

3. Acceleration effects can be reduced by reducing the particle size 
of the. metal additive, the oxidizer particle size, and the metal 
content. 

4. The effect of acceleration on burning rate is given by the equation 


Ib_ r l-n(o-i) 1 cos * 

(r ) L 1 - II [1 + {$ ~ coa $ - 1) G ] J 

' b a=o $ = o 

where G is a function of the variables ( 7 and (ap /r^) /(a p / r^) c . 

5. The theoretical model is in qualitative agreement with experiment 
with regards to the effect of acceleration on the pressure exponent 
and the burning rate coefficient and the effect of non-normal 
acceleration forces on burning rate. 

Nonmetallized Propellants 


1. The pressure exponent is not strongly dependent upon acceleration. 

2. There is an upper bound on acceleration induced changes in burn- 
ing rate. 
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3. Burning rate is insensitive to accelerations under 100 g. 

4, The granular diffusion flame model predicts the following: 


a. An upper bound on burning rate. 

b. The burning rate ratio is independent of initial tempera- 
ture and catalyst content when the acceleration force is 
parallel to the burning surface, 

c. Acceleration effects can be reduced by adding burning 
rate catalysts, reducing the mean size of the oxidizer 
particles, and increasing the initial temperature of the 
propellant. 

d. The burning rate ratio is given by the equation 



a=o 


Gr 


1/2 


cos # 


~l/2_ 

C _ Re 
d o 



Gr cos $ 
d 


C Re 2 
d o 



o 


5, The effects of acceleration induced pressure differences across 
and mean flows in the gas phase reaction zone are negligible at 
accelerations under 1000 g. 


6. The major effect of acceleration is derived from the effect of 
acceleration on the heterogeneous structure of the gas phase 
reaction zone. 


7. The theory correlates the experimental data of Reference 37 for 
accelerations above 200 g. 



RECOMMENDATIONS 


The following recommendations are made, 

1. Internal gas dynamics studies of motors with slotted tube grains 
(circumferential slots) and CP grains with uninhibited ends should 
be made because large swirl velocities can be achieved from 
radius changes; these swirl velocities can affect burning rate 
(erosive burning) and metal/metal oxide retention. 

2. A theory for particle burning based solely on heat transfer (no 
metal combustion) should be developed, 

3. The theory for metallized propellants developed herein should be 
examined to determine the effects of acceleration induced 0 
changes. 

4. A more general combustion model for the heterogeneous gas phase 
reaction zone should be formulated, 

5. Experimental data for the thrust and discharge characteristics of 
nozzles with swirling flows should be obtained. 

6. The effects of non-uniform burning rates on pressure and thrust- 
time histories should be determined, 

7. Additional data concerning the effects of acceleration on the burning 
rate of both metallized and nonmetallized composite propellants 
should be obtained. 
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APPENDIX A 

EXPRESSIONS FOR dip / ip AND d <p /<p 


Consider the function 
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where a is a constant. It is readily recognized that the function g has the same 
form as Ip and tp depending on the value of a . Application of Liebnitz Rule yields 
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APPENDIX B 
FLOW IN A CONICAL PIT 


The steady, axis ymme trie, and inviscid flow of an incompressible fluid in a coni- 
cal pit was studied. Figure B-l is a schematic diagram that illustrates the cone, 
coordinate system, and nomenclature. In addition to the assumptions mentioned above, 
it was also assumed that the combustion gases leave the surface of the cone along a 
normal at the velocity 


V e (r, 8.) 


r, p /p 
b P 
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and that the cone is of infinite extent. 

The axisymmetric assumption is based on the symmetry of the problem and the 
compressible assumption is based on the fact that the Mach numbers will be small 
since the acoustic speed is large. The assumption concerning the boundary condition 
is based on photographs of the flow in cones around fibers. The inviscid and infinite 
extent assumptions are made to simplify the problem. 

With the assumptions above, the governing equations are 
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and the boundary conditions are Eq, (B-l) and 


V Q (r, o) = V (r, 0 ) = 0 (B-5) 

a r c 

The pressure terms in Eqs. (B -3) and (B-4) were eliminated by cross differentia- 
tion and then summation. The resulting equation and Eq. (B-2) were reduced to ordinary 
differential equations by assuming that the dependent variables were functions of theta 
alone . 


This yielded 
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2 f + g + g cot 9 = 0 


gf + f g = = g g = 0 


The boundary conditions became 


f (0 c ) = g (0) = 0 
g( 0 c } = - Vp /p g 
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A numerical technique employing the Runge-Kutta method to integrate from 0 = 

6 to 0 = Aand a Taylor's series to extend the results to 0 = 0 was developed. The 
resulting computer program was partially debugged. 
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APPENDIX C 

TYPICAL PROPERTIES OF NONMETALLIZED 
COMPOSITE PROPELLANTS 


An important portion of the analysis presented herein is concerned with the 
estimation of the order of magnitude of various effects. For these estimates, typical 
values for various thermo-physical properties are required. The objective of this 
appendix is to present typical values for the required properties. 

Nonmetallized composite propellants essentially consist only of binder and 
oxidizer. The binder is usually based on a hydrocarbon material; the pertinent 
oxidizer here is ammonium perchlorate. The general characteristics of the propellant 
are selected in practice as a balance between energetic and physical property considera- 
tions. In most cases of interest, the weight fraction of binder ranges from 10 to 25 
percent. 

Figure C-l presents the effect of the weight fraction of binder on the adiabatic 
flame temperature, molecular weight, and specific heat ratio of the products of com- 
bustion of the propellant at equilibrium at a pressure of 1000 psia. The results were 
calculated by the method of Reference 43 and employed JANAF thermophysical 
properties. The figures shows that typical values for these variables are as follows: 

T f = 4500 R, y = 1. 22, and M = 22. 

Reference 44 has measured the surface temperature of nonmetallized composite 
propellants with ammonium perchlorate oxidizer and has found that the surface tempera- 
ture ranges from 1500 to 1700 R. Thus, a typical value for the temperature of the 
burning surface is 1500 R. 

The mean temperature of the gas in the gas phase reaction zone is the mean of 
the surface and flame temperatures. Therefore, a typical value for this temperature 
is 3000 R. 

Figure C-2 presents the effect of temperature on the thermodynamic properties 
of the pyrolysis products of a typical binder at 1000 psia. The binder is the one 
employed in the previous calculations and the computational procedure employed is 
the same. The figure shows that typical values for the molecular weight and specific 
heat ratio at the mean temperature of the gas phase reaction zone are 19 and 1. 12, 
respectively. 

Figure C -3 presents the effect of temperature on the thermodynamic properties 
of the equilibrium decomposition products of ammonium perchlorate at 1000 psia. The 
computation procedure is the one employed previously. The figure shows that the 
molecular weight and specific heat ratio at the mean temperature of the gas phase 
reaction zone are 28 and 1.24, respectively. 

It is important to note that the molecular weight and specific heat ratio of the 
fuel and oxidizer vapors are not widely different from those of the products of combus- 
tion. Therefore, these properties are relatively constant throughout the gas phase 
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reaction zone. The values given previously for the products of combustion of the pro- 
pellant will be employed as typical values. 

Typical values for viscosity, Prandtl number, and thermal conductivity were 
computed from the previous typical values and the relationships presented by Bartz^. 
Thus, fig = 1 x 10“6 slug/ft-sec and Pr^ = 0.82, 

The densities required were computed from the perfect gas law for the previous 
typical conditions (p = 1000 psia, M = 22, Tf = 4500 R, t s = 1500 R, t g = 3000 R). 

This yields f> 3 = 0 . 042 slug/cu. ft. , P = 0 . 021 slug/cu. ft. , and P £ = 0 . 014 slug/ 
cu. ft. 

The burning rate of nonmetallized composite propellants usually ranges from 
0. 2 to 2.0 in/ sec. Therefore, a value of 1 in/sec was chosen as typical. 

The density of nonaluminized composite propellants is around 0.060 lb m /cu. in. 
This value will be employed herein. 

The normal velocity of the gas at the burning surface can be computed from the 
continuity equation or 


v s = VV P S 

Substitution of the pertinent typical values gives 

v g = 6.4 ft /sec 

Reference 32 has employed several different methods to measure the thickness 
of the gas phase reaction zone of a nonmetallized composite propellant with ammonium 
perchlorate oxidizer. All of these methods indicate that the zone is less than 100 
microns thick at pressures above 200 psia. Therefore, it will be assumed that 100 
microns is a typical value for the thickness of the gas phase reaction zone. 

A typical value for the mean diameter of the oxidizer particles is 50 microns. 

The typical values computed and assumed are presented in Table C -1. 
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TABLE C-I 

TYPICAL VALUES FOR THERMOPHYSICAL PROPERTIES 


p 

= 

1000 psia 


= 

1 x 10“ 6 slug/ft-sec 

T f 

= 

4500 R 

P s 

= 

4. 2 x 10-2 slug/ft^ 


= 

1500 R 


= 

2. 1 x 10"2 slug/ft^ 


= 

3000 R 

P f 

= 

1.4 x 10-2 slug/ft^ 

M 

= 

22 

P P 

= 

3. 2 slug /ft ^ 

v 

= 

1. 22 

v s 

= 

6. 4 ft /sec 

Pr 

g 

= 

0. 82 

v f 


19. 1 ft /sec 

R 

= 

2.27 x 10*^ ft -lbf /slug -R 

r b 

= 

8.33 x 1 0“ 2 ft/sec 

^ox 

= 

1. 64 x 10* 4 ft 

6 

r 

= 

3. 28 x 10' 4 ft 
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